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ABS TRACT 


Improvement of water supply and water quality in the Cooking 
Lake Moraine near Edmonton is the key to more intensive recreational 
use and other uses of the area. One of the ways of improving both might 
be found in better snow management. To accomplish this an improved 
method for estimating regional and local patterns of snowfall is needed 
so that snowmelt runoff can be estimated. 

There are three main parts in the study. The first is a review 
of the landform, vegetative cover, meteorological data, and other factors 
that might affect snowfall and snow drifting patterns. The second 
involves the use of LANDSAT I imagery and ground truthing in various ways 
to map snowfall and snow drifting patterns on the ascent: The third 
part contains an analysis of these patterns using Thornthwaite techniques 
in determining local and regional water balance patterns. 

Snow depth and density studies were conducted at 12 sites with 
a Mt. Rose snow sampler. These sites were identified on photographic 
plates. 

Based on the information obtained from the snow surveys, 
Thornthwaite water balance calculations, and LANDSAT I imagery it was 
possible to estimate runoff potential. Calculations of potential snowmelt 


runoff were made for several areas on the moraine. 


1 ANDSAT I imagery was discussed under five subheadings: Scale 
1:£,000,000 LANDSAT I Imagery, Scale 1:250,000 LANDSAT I Imagery, 
Density Sliced LANDSAT I Imagery, Tri-Color Projected LANDSAT I Imagery 
and Water Balance Relationships and LANDSAT I Imagery. 
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The procedures employed in this thesis, with the modifications 
suggested, should be useful in future investigations of snow distribution 
patterns. With more in-depth information about these patterns it may be 
possible to manage snowmelt runoff for both water supply and water 
quality in the moraine and in other areas with similar problems and 


Opportunities. 
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Chapter I 


INTRODUCTION 


One of the most dynamic geophysical features of Planet 
Earth is the ever-changing snow-cover. No other wide- 
spread soil material on the earth's surface is deposited-- 
and is removed--so quickly (Meier, 1973, p. 1). 

The purpose of this study was to develop a more accurate means 
of estimating regional patterns of snowmelt runoff using remote sensing 
techniques. The area selected for study was the Cooking Lake Moraine 
twenty-five miles east and south-east of Edmonton (Plate 1). The 
area is a hummocky dead-ice moraine with widespread sloughs, minor 
depressions and deep forest soils. 

There are major water quantity and quality problems relative 
to recreational facility development and to other uses of the Cooking 
Lake Moraine. Restoration of the lakes on the moraine to levels and 
qualities experienced in a relatively wet period early in this century 
is a widely supported objective. A pipeline from the North Saskatchewan 
River (near Devon) to Miquelon Lakes and a cascading system to other 
lakes was one of the proposed solutions to the water aan piiey problems 
(EPEC, 1970). Pipelines could be built and maintained at a substantial 
cost to the public. A less expensive alternative might be found in 


better management of regional snowmelt runoff each spring and in use 


of regional or nearby surpluses. 
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Legend for Plate 1 


A December 5, 1974 1:1,000,000 LANDSAT I 


Image of the Edmonton Area 
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The objectives for chee study were: 

1. To establish what variations in winter snowfall are present in the 
Cooking Lake Moraine area; 

2. To develop a procedure using LANDSAT I imagery for determining 
snow concentration and distribution patterns; 

3. To investigate the relationship between LANDSAT I images' tonal 
variations and snow concentration and distribution patterns; 

4, To experiment with various enhancement devices to maximize dis- 
crimination of snow cover variations; 

5. To compare regional snow moisture recharge capacities and snow 
water content levels required to produce runoff; 

6. To investigate the relationship of snow concentration and distri- 
bution patterns to soil moisture storage capacities so that runoff 
potential might be estimated. 

Interpretations from the LANDSAT I imagery were made with two 
major assumptions in mind. The first assumption was that snow concen- 
tration and distribution patterns would not vary appreciably from 
January 10, 1975 to April 10, 1975. This was a necessary assumption 
because the April 10, 1975 LANDSAT I imagery could not be obtained 
(Refer to Chapter V for details). The second assumption was that 
meteorological data for the Edmonton International Airport combined 
with 1975 winter snow depth and water content measurements taken with a 
Mt. Rose snow sampler would be representative of the major patterns 
in the Cooking Lake Moraine. 


The bases for water quality and quantity studies for the 
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Cooking Lake Moraine developed out of public concern for lake improve- 
ment. Historically the Cooking Lake Moraine has experienced rather 
large fluctuations in water quality and quantity. In some instances 
man has played a major role in these fluctuations. The following 
chapter is a review of studies elsewhere that might serve as bases 


for the analysis of these problems. 
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Chapter II 


REVIEW OF LITERATURE 


This chapter is a review of the literature pertaining to the 
Cooking Lake Moraine, and an explanation of the enhancement techniques 
used to show snow concentration and distribution patterns with LANDSAT I 
imagery. The chapter is subdivided into four parts: Geological 
History, Settlement, Lake Levels, and LANDSAT I Imagery and Snow 


Drifting. 


Geological History 


"During the Pleistocene epoch, a continental ice sheet 
originating on the Precambrian Shield in Keewatin Territory advanced 
from the northeast and covered the region at least aetage are 
1962, marginal notes on Surficial Geology map 83H). The second 
advance and recession, in late Wisonsin time, formed most of the local 
topographic features of the Cooking Lake Moraine. The Cooking Lake 
Moraine is a "hummocky moraine: .. . till composed of mixed clay, 
silt and sand, with pebbles and boulders; lenses of sand, gravel and 
local bedrock; generally more than 40 feet thick; topography undulating 
to gently rolling" (Bayrock, 1962, marginal notes on Surfical Geology 
map 83H, Figure 1). The area to the south-west of Joseph Lake, is a 


"oround moraine: .. . till composed of mixed clay, silt and sand, 
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with pebbles and boulders, variable in thickness, but generally less 
than 40 feet, topography level to undulating" (Bayrock, 1962, marginal 
notes on Surficial Geology map 83H, Figure 1). Most of the Edmonton 
district surficial materials are formed on a bedrock surface sculptured 
in poorly consolidated materials, which comprise the Edmonton Formation 
of late Cretaceous age (Carlson, 1967, p.-5). An erosional remnant of 
this material served as a base for the Cooking Lake Moraine and provided 
much of the regional relief. The district was comprised of three main 
formations: the Edmonton, the Bearpaw, and the Pale and Variegated 
(Figure II). Carlson (1969) provided an in-depth analysis of the 
bedrock topography in ''Bedrock Topography and Surficial Aquifers of 
the Edmonton District, Alberta''. Since glaciation there have been very 
few topographical changes on the moraine. The relatively dry climate 
combined with the lack of a major river has helped to keep erosion and 
deposition to a minimum. Minor filling, though, has occurred with 
vegetation, in most of the sloughs. 

Cooking Lake loam is the predominant soil of the area 
(Figure III). This loam and other soils such as Uncas loam, Falun 
loam, Angus Ridge loam, Camrose loam, and others are all present on 
the Cooking Lake moraine. For an in-depth study of these soils refer 
to the Soil Survey of Edmonton Sheet (83H) 1962, by Bowser, 
Kjearsgarrd, Peters, and Wells; Grey Wooded Soils and Their Management 
by the University of Alberta, Faculty of Agriculture; and Canadian 


Soil Science by Atkinson. 
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Source: Surfical Geology Map Edmonton N.T.S. 83H 1972 
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Source: Soil Survey of the Edmonton Sheet 83H 
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Settlement 


"About 18 miles east of Edmonton lies a hilly area which some 
people call the Cooking ‘ake Moraine, while others use the old Indian 
name Beaver Hills (Amisk Wachee)" (Nyland, 1969, p. 22). It was 
assumed that the origin of the name "Beaver Hills" came from the Beaver 
Indians who inhabited the area. About 1780, the Blackfeet Indians 
passed into the Beaver Hills region during their migration from the 
south eastern prairies, and later the Crees invaded and forced the 
Blackfeet Indians from the Beaver Hills. There was an abundance of 
game--bison, moose, elk, deer, ducks, geese, and swans--on the Cooking 
Lake Moraine. Those rich hunting grounds led to a number of conflicts 
among the various tribes living in the area. 

By 1895, the Cooking Lake Moraine was sparsely inhabited by 
white settlers. They cleared the land and burned the trees. Some of 
the fires were blown out of control, destroying large stands of natural 
forest. The Edmonton Bulletin described the destruction on May 9, 

LOO J 
The worst of the matter is that many settlers look upon this 
wholesale destruction of timber as something commendable and 
to the advantage of the country in making more land available 
for settlement. Even in those parts of the country that are 
open, only a very small percentage of the land is actually 
cultivated, so that there is no necessity for destroying the 
timber wholesale. (Nyland, 1969, p. 22) 
According to Nyland the rain had put the fire out in August of the same 
year. In October of 1895, the fire was reignited and the flames 


crossed the North Saskatchewan River, burning on northward; leaving 


only small stands of trees. Some of the settlers had a more constructive 
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idea for the natural surroundings. The first tree replanting conser- 

vation practices originated in 1892; by 1898 the first ranger station 

on the moraine was established and during 1910, the first tree nursery 
in Alberta was established on the shores of Cooking Lake (Nyland, 


1969) 


Lake Levels 


"In 1865 both Beaverhill Lake and the lakes in the moraine 
were dry, but heavy rains in 1875 restored the water levels'' (Environment 
Conservation Authority, 1971, p. 2). Appendix A, indicates the 
progression of lake level fluctuations as calculated by Laycock in an 
article entitled "Lake Level: Fluctuation and Climatic Variation in 
Central Alberta" (1973). The years 1883 and 1888 were the only years 
prior to 1900 which showed a surplus and a positive lake level change. 
There are conflicts in the literature concerning lake level fluctuations 
from 1904-1930. Research by Acres West Limited indicated that the 
lakes were in satisfactory condition for most uses until the 1930's 
when water teres declined steadily and the lake quality began to 
deteriorate (Acres Western Limited, 1971). The surplus and 
deficiency patterns presented by Laycock, Appendix A, showed that 
from 1900 to 1904, there was a major rise in lake levels. Also, 
surpluses; inei907,, 1914, 1917, 1920, 1925, and 1927 combined with a 
negative or very minor lake level change during the remaining nineteen 
years (1900-1930) indicated that neither a stationary pattern nor 


inadequate lake levels existed for most recreational activities. 
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Stanley Associates Engineering Limited in a report entitled "Water 
Inventory and Demand Draft Report", (1974) agreed with Acres Western 
Limited, that the lakes were in good condition before 1930. 

It was apparent from reviewing EPEC's and Stanley Associate's 
lake level graphs that the lakes on the Cooking Lake Moraine have been 
fluctuating continuously. Both EPEC and Stanley Associates used these 
graphs to indicate a continual recession in lake levels, even though 
data used to construct these graphs were very limited. For example 
most of the data for Cooking Lake were estimated from interviews with 
local residents. For Ministik Lakes less than two years of data were 
available thus rendering any conclusions made from these graphs to 
describe lake patterns in the 1930's meaningless. 

For Miquelon Lake, negative lake level fluctuations have 
occurred as a direct result of man's activities. 

No mention of a noticeable drop in lake levels was made until 
1925. A settler had slightly deepened the outlet of the lake, 
but the big blow did not come until 1927 when a water shortage 
in the Camrose town reservoir was temporarily solved by digging 
a canal through the height of land south of Miquelon Lake to 
divert water into Stony Creek. There appears to have been no 
permission given for this diversion. The move was, to say the 
least. Hu advisede(Nyland; 1970, epin22)is 
In reality the small drainage basin of the Miquelon Lakes suffered 
several dry years (approximately 10) and thus was unable to withstand 
the drainage draft put on it by the ditch into the Camrose Creek Basin. 
"When in 1928 the lake level dropped lower than the control weir, the 
canal was dug deeper with the help of heavy machinery until on the 


height of the land the canal was close to 20 feet deep" (Nyland, 1970, 


p- 22). In 1971, Nyland went on to say, that Calgary Power had dug’a 
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ditch which drained Miquelon Lakes 16 feet (Nyland, 1971, Public 
Hearings). Calgary Power could not be blamed for the entire 16 foot 
drop in lake levels, particularly when lake levels dropped well below 
the ditch during the next few years. 
Deforestation caused by fires and land clearing has had an 
effect on land use on the Cooking Lake Moraine. In many areas the 
"parkland" was changed into grain fields and pasture land. The effects 
of deforestation on water retention capacities have been debated. 
Deforestation and topographical changes associated with land 
cultivation and the advent of modern civilization in the water- 
shed and surrounding areas has resulted in the deterioration of 
water retention capabilities of the watershed, consequently 
ground water storage and water tables in the watershed have been 
Patino (Zuzak o/h, Public Hearings, | p. 7). 

Zuzak also accepted Nyland's explanations about falling lake levels. 

Nyland named land clearing operations and slough drainage projects 

in the watershed as having adverse effects on runoff. "These previously 

timbered lands used to provide a reliable source of water for the 

springs in the lake, which now flow only intermittently" (Nyland, 

1970, p. 25). In the author's opinion, deforestation would add to 

lake levels due to a decrease in soil moisture storage utilization 

by plants. Drops in yearly precipitation patterns, in this case, 

appear to have had a greater effect on lake levels than did deforestation. 
It has been suggested that the clearing and burning of forest 
cover in the Cooking Lake basin has resulted in lower flow into 
the lake. In principle, such clearing and burning would result 
in a reduced soil moisture consumption and thus in yield 
increases, particularly in a sub-humid region such as this. In 
practice, the extensive burning before 1900 probably contributed 


to the very heavy yields of the wet years that followed because 
of the reduction in consumptive use (Laycock, 1973, p. 89). 
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Similar variations in thought (non professional versus professional ) 
typified much of the discussion in the literature available about 
the Cooking Lake Moraine. 

Major interest in the water quality and quantity problems 
of the Cooking Lake Moraine resulted in a series of studies beginning 
in the late 1960's. In 1969, the Conservation and Utilization Committee 
of the Alberta Department of the Environment in cooperation with the 
Cooking Lake Area Study Committee established a conservation area. 

In the fall of the same year Nyland presented a description of the 
Moraine's settlement, fires, forestry, recreation, precipitation, 
and water levels in an article entitled "The Dying Watershed". The 
article provided a detailed review of the history of the Cooking Lake 
Moraine, but lacked reference to in-depth water studies which would be 
required for the type of conclusions which were drawn. 
The watershed can be restored, then the beaches will come to 
life again and the fish will return to be enjoyed by man and 
beast. Grazing areas can be placed under proper management 
and possibly support more cattle than is now the case. Those 
soils most suitable to trees can be replanted to the benefit 
of wildlife, for recreation and to help restore the balance 
of nature. Both the people of the city and the rural residents 
can only benefit by such a scheme (Nyland, 1969, p. 38). 

During the remainder of 1969 and 1970 a number of reports were 
published concerning the Cooking Lake Moraine. Many of the publications 
dealing with the Cooking Lake Moraine were repetitive and would add 
little to this review. 

Environmental Planning and Engineering Consultants (EPEC) 


of Edmonton, Alberta produced the first consultant's report on June 10, 


1971, entitled "Economic Analysis of the Cooking-Hastings Lakes". 
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The report was developed at the request of the Water Resources Division, 

Department of Agriculture, Province of Alberta. 
More specifically the study was to incorporate the investigation 
of three areas. Firstly, a verification of the authenticity 
of a petition calling for "Action to Reclaim the Beaverhill 
Watershed" presented to the Environmental Conservation Authority 
in ithnestallyot 29/70. “This was. a requirement, in: order ‘to be able 
to determine the validity of the petitioner's alleged concern 
respecting Cooking and Hasting Lakes. Secondly, an investigation 
as to the possibilities of improving the conditions of the lakes 
through internal water management and/or the importation of water 
from sources outside the Cooking-Hastings watershed. This 
necessitated an examination of possible schemes and their rami- 
fications. Thirdly, an economic analysis identifying probable 
costs of possible schemes and identifying benefits that could 
accrue in various sectors such as recreation as a result of 
improvement to the existing conditions of the lakes (EPEC, 1971, 
Pe l-la). 

During the fall of 1970, a total of 598 people signed a 
petition calling for "Action to Reclaim the Beaverhill Watershed" (EPEC, 
1971, p. 5). Only 284 signatures were available and of those only 222 
had legible names and addresses. A questionaire was sent to 25 
individuals on the list: 24 people completed the questionaire. The 
results of this questionaire were presented to EPEC by Bailey, Director 
of the Water Resources Division, (Appendix A page vii of the EPEC report). 
Unfortunately or fortunately depending on how the results were to be 
used it was found that twenty of those petitioned owned property in the 
Cooking Lake Moraine. Most property owners would enjoy an opportunity 
to increase their property values. The prices were increasing and will 
undoubtedly continue to grow as development occurs (Table 1). In all 


instances the mean increase in purchase price was 36.2 percent 


and the properties were resold in an average of 13.5 months. 
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Table 1 


PRICE INCREASES FOR SELECTED PROPERTIES 


IN STRATHCONA COUNTY 


Percentage 
Purchase Selling Holding Increase Over 
Location Size Price Brice Period Purchase Price 
(Tp./Range ) (Acres) (s) (Ss) (Months) 
SLIZO 161 8,700 16, 000 38 84 
sal iealt 161 12,000 18 ,000 38 50 
SLI ZL 161 85390 10,000 23 ib 
oh hea 161 2.500 3, 668 1 47 
sy hy Wa 161 lle Sia By 133/00 6 (e: 
Bi a: 166 10,000 16,000 19 60 
sid Ya 601 67,000 87,500 14 31 
SL i22 40 30,000 35,000 3 i 
Sy ay are 6 Syysle\s) 77000 1 36 
a 3 4,400 6,500 3} 48 
ey Pas 38 11,000 12,000 Zz 9 
53/22 20 10,000 10,800 8 8 
payee 20 6,000 6,800 Zee V3 
54/20 L59 55900 9,500 29 13 
54/20 161 3, 700 5,000 1S) 35 


Source: Acres Western Limited, 1971 "Alternative Land Use Evaluation 
Cooking Lake and Hastings Lake Area". 
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In many instances EPEC seemed to lack adequate hydrological 
information. 'Water balance relationships are much better established 
and the effects of the watershed management are much better known pro- 
fessionally, than is suggested in the EPEC report" (Laycock, 1971, 
Public Hearings, p. 79). 

Hydrological matters were not the only subjects in the EPEC 
report to receive criticism. "We are of the opinion that the economic 
analysis of Cooking and Hastings Lakes prepared by EPEC Consulting 
Ltd. should not be the bases for the decision on whether or not to 
invest in restoring the water levels of these lakes" (Lefrancois, 
UA gm Cat No Bien eariae ane. 46). The report was narrowly mltiple- 
purpose with consideration of recreation, cottage development, residential 
development, commercial devetopment, and waterfowl needs. EPEC's 
Strategy, which aimed for multiple means, but did not achieve them for 
raising lake levels, was represented by three pipeline schemes (A, A', 
and B) for pumping water from the North Saskatchewan River to various 
points on the Cooking Lake Moraine (Figure IV). In 1971, the direct 
costs were estimated to be $4,128,840 for scheme A, $4,408,680 for 
scheme A', and $7,971,160 for scheme B. If a choice of one of three 
pipeline schemes had been made in 1971, it appeared from a review 
of the Public Hearings that scheme B would have been favored. If this 
decision had been made it would have been based almost totally on 
estimations rather than fact. Mayon Dent of Edmonton stated, "In our 
opinion, the recreational benefits would be dramatically increased 


by using scheme B over scheme A'' (Mayor Dent, 1971, Public Hearings, 
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p- 13). The Fish and Game Department endorsed scheme B (Chymko, 1971, 
Public Hearings). "Tentatively however, I would opt for No. 6 

in, the tact that (have, or plan’ "B’in the BPEC report™ (Laycock, 

1971, Public Hearings, p. 83). The place of intake for scheme B near 
Devon would result in fewer pollutants being present than for schemes 

A or A'. "The problem of high bacteria counts in the North Saskatchewan 
River downstream of Edmonton is well known, and it has been the subject 
of much discussion already at these Hearings" (Simpson, 1971, Public 
Hearings, p. 94). 

Four years later Gallup (1975), after conducting intensive 
water quality studies, stated that in raising the level of water on 
the Cooking Lake Moraine the only effect would be to create large 
sloughs out of small ones (Gallup, 1975, Cooking Lake Management 
Committee Meeting). In the same year Underwood, McLellan and Associates 
stated, "That while the level of Cooking Lake may be raised, giving 
more area for boating, the water quality won't be changed much"! 
(Edmonton Journal, 1975) 

Laycock, in 1971, and in 1973, suggested internal modification 
of lake levels using drainage of one of the larger lakes such as Joseph 
Lake and of local depressions as a means of restoring Cooking Lake 
and Hastings Lake. Laycock indicated that the cost for this type of 
internal management would be approximately one tenth of that of the 
pipeline schemes. In a 1973 article, several other alternatives were 
suggested by Laycock, one of these was to pump water from adjoining 


basins, such as Little and Big Hay Lakes, in the spring into the Cooking 
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Lake Basin and using Looking Back Lake and other water bodies as storage 
reservoirs to extend pumping seasons. 

In light of Gallup's and Underwood, McLellan and Associate's 
statements, spending approximately $8 million for massive pumping from 
the North Saskatchewan River did not appear to be the best alternative. 
Also, if construction cost increases parallel many of the increases 
for the Corp of Engineers projects in the United States, or dredging 
Operations on the St. Lawrence River, the costs may double or triple 
before completion of the project. 

The second consultants report was presented to the Cooking 
Lake Area Study Management Committee on October 1, 1974, by Stanley 
Associates Engineering Ltd. The report was entitled "Water Inventory 
and Demand--Draft Report''. The Stanley Associates' Report came closer 
to developing a multiple purpose and multiple means approach to the water 
quality and quantity problems which the Cooking Lake Moraine now endures. 
The Stanley Associates study included research, but there are many 
instances where they employed research as a tool without adequate thought. 
"A detailed analysis was carried out on Hastings, Cooking, Ministik, 
and Miquelon Lakes to determine the correlation, if any, between rain- 
fall and lake levels for 1972" (Stanley Associates Engineering Ltd., 
1974, p. 40). Three years earlier at the Public Hearings, Kellerhals 
said that five years of hydrological data are needed before conclusions 
could be drawn (Kellerhals, 1971, Public Hearings). The Stanley 
Associates report was, in many instances, based on only one year's 


data. Two internal drainage schemes were presented by the Stanley 
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Associates, the draining of sloughs and other small pockets of water 
and the draining of one of the lakes in the moraine. Their conclusions 
were, "In view of the relatively large water volume deficit calculated 
for the moraine it is probable that such a scheme would only provide 

a small percentage of the volume of water required to just stabilize 
the lakes at their present level, and have a high cost benefit ratio" 
(Stanley Associates Engineering Ltd., 1974, p. 103). In some years 
this may be true, but in years of above average precipitation, such as 
1974, large volumes of snowmelt runoff would be available for internal 
management schemes. 

The cost of internal water management would be small compared 
to piping water from the North Saskatchewan River. In most years 
stablization of lake levels (for Cooking and Hastings Lakes) could be 
obtained by internal management of snowmelt runoff and by draining 
some sloughs and small lakes. 

The purpose of the report was centered around the Dodds 
Thermal Plant. 

It has not been possible to carry out the in-depth study we 
originally envisaged. This was due to two reasons. Firstly, 
time constraints were imposed on the study by the need to 
investigate the possibility of a multi-purpose pipeline being 
constructed to serve the needs of the proposed Dodds Thermal 
Plant and the Cooking Lake Moraine. The time constraints 
{were} being imposed by the need for Calgary Power to have 

a decision on such. a scheme by March 1, 1975. Secondly, 
paucity of data has precluded in-depth analysis of the 
situation with respect to the hydrology of the Cooking Lake 
Moraine (Stanley Associates, 1974, p. 2). 


Scheme B from the original EPEC report had been altered to include 


the Dodds Thermal Plant. "One gets the impression that both may have 
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been influenced by a prior knowledge of government desires to have a 
strong case for a pipeline through the area to serve thermal power 
development to the east" (Laycock, 1974, p. 10). 

At a meeting of the Study Management Committee, October ll, 
1974, Laycock, the co-ordinator of the University of Alberta's Water 
Resource Center Lake Area Study Program, presented objections pointing 
out the inadequacies in the report. Deeprose, the Chairman of the 
meeting, suggested that Laycock be contacted concerning elaboration of 
points raised at the meeting. He was contacted once, a week later, 
concerning groundwater aspects of the report. Laycock then prepared 
a paper highlighting certain aspects of the report and presented it to 
Stanley Associates. "Our Study Program group was not satisfied with 
this degree of contact, at least partly because we had supplied some 
of the background research data (almost unacknowledged) which we 
believed had been ill-used in the interpretations and conclusions 
piven mm the dratt report” (Laycock, 1974) p. Pry. 

If lake levels were to be raised, the cheapest, most efficient 
methed should be employed. In the author's opinion better management 
of snowmelt runoff and the sacrificing of one or more of the lesser 
used lakes would provide a viable solution in most years. If more 
water was needed the Little Hay Lake drain flow to the south of the 
Cooking Lake Moraine, could be used in the spring of each year. This 
could be complemented with flow from the Clearwater Creek and the Irvine 
Creek Basin from the west, and the Katchemut Creek Basin from the east, 
which could be pumped to various points on the moraine. Thus the 


purpose of this thesis was to develop a more accurate means of 
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estimating regional patterns of snowmelt runoff using remote sensing 
on and around the Cooking Lake Moraine. LANDSAT I imagery proved to 


be an efficient tool from which snowmelt runoff could be estimated. 


LANDSAT I Imagery and Snow Drifting 


"Snow cover is the primary source of Canadian water supplies" 
(McKay and Findlay, 1971, p. 17). Snowmelt and spring rains are 
SUbstantialsources Of gunotr (Laycock, 19/1, Public Hearings, ps 80). 
The amount of snowmelt runoff, in most years, determines much of lake 
level fluctuations on the Cooking Lake Moraine. Obviously there are 
other factors affecting lake level fluctuations, such as summer rains 
and evaporation, but in average and above average snowfall years snow- 
melt runoff causes the greatest lake level fluctuation. 

To properly manage this resource, there is a need for 
continual synoptic data on snow distribution patterns. "Continuous 
monitoring of the highly variable snowcover is possible only from 
satellites" (Meier, 1973, p. 1). Research by Wendler and Carlson (1974) 
supported the contention that LANDSAT I imagery could be used success- 
fully to monitor snowmelt. 

The major complication with LANDSAT I imagery Leas cloud cover. 
Yet, there are a number of images available for the fall and winter 
periods of each year except 1974, when only two cloud free images were 
available for the Cooking Lake Moraine, one in November of 1973 and one 


in April of 1974. 


ria bats ee 
ang Jee ray ‘“, 
7 sn Wig 


i Ws tea ia 
pat 4 q 
as x) ert » * ei 4 5 —— 
‘ Ps r ' AVR She, ms i 
rie Dey Crue ge ; i atyar's © 
eT) ‘ 
j 
! ' 
i 
arn a A: 
i 


i 
re ” 
- Atha 
| { LaF RS ee. 
7 i : : 
: i vad 
itihiy bid a4 — ” 
= ‘ : i? - % Fe a . ’ 
‘ ay 
sare Seth 


se 
ard Nenad a) 
& Uy 


A number of snow LE have been conducted using remote 
sensing techniques. Meiers (1973) evaluated LANDSAT I imagery for 
mapping and detection of changes in snow cover on iand and on glaciers 
(Northern Cascades in Washington). Meiers found that under favorable 
conditions the mapping procedure used was reliable to within four 
percent of the snow covered area. Meiers (1973) monitored drainage 
basins by using the Stanford Research Institute console to electronically 
superimpose basin outlines on imagery and by using video density slicing 
to measure areas. A zoom transfer scope was also used to transfer 
information from 1:250,000 prints to standard topographical maps. 

Barnes (1973) evaluated the applications of LANDSAT I imagery 
for mapping snow cover in mountainous regions of the western United 
States (Salt River Project in Arizona, 1973). Barnes found that band 
5 (.6 to .7 um) was the most useful band for detecting and mapping 
mountain snow cover (Barnes, 1973). Barnes also found that the 
near infrared band (7) proved to be useful in detecting areas of 
melting snow. 

The hydrologic characteristics of snow covered terrain with 
thermal infrared imagery (Near Cape Discovery, Ellesmere Island) were 
investigated by Poulin (1973). "The energy that penetrates the snow 
over land is either absorbed or reflected at the snow/soil interface, 
while at a snow/ice interface a significant portion of energy continues 
to be transmitted into and through the ice" (Poulin, 1973, p. 5). 
Several graphs were presented showing the thermal conductivity of snow 
compared to sea ice, fresh ice, frozen silt and clay, and frozen 


sandy silt to add significance to his previous statement. Poulin 
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indicated that a satellite would be the preferred vehicle for obtaining 
thermal imagery, but with present systems the noise level would cause 
large temperature differentials. 

The majority of the articles published on remote sensing and 
snow surveys deal with the mountain environment; consequently many of 
the investigators found problems with slope angles, cloud cover, snow 
depth, and accessibility. A complete list of publications with 
references to the use of remote sensing for snow surveying was provided. 
by the Canada Centre for Remote Sensing (Appendix B). 

Based on the review of literature, several remote sensing 
systems were found to be superior for snow surveys. A radar system 
and a microwave radiometer system (antenna systems) produced the most 
useful images and printouts available for snow surveys. These systems 
could be used to display the general snow cover and could be adjusted 
to penetrate the same snow cover to calculate snow depth (water content 
measurements would have to be made in the field). A thermal infrared 
scanner (optical system) could be used to provide detailed information 
about varying degrees of snow ripeness and water content. Computer 
printouts from LANDSAT magnetic tapes provided a digitized format of 
the study area. The "Imager 100'' sorts digitized data and is capable 
of identifying the spectral signatures of a target, and of sorting 
the digital data to discriminate and provide an image prediction, through 
correlation analysis, for the occurrence of a particular target. 

The remote sensing systems listed represent some of the finest 


presently available, in capabilities for a snow survey and also some of 
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the most expensive to use.The microwave radiometer, radar, and thermal 
infrared systems are presently used only on aircraft platforms, making 
continual coverage very expensive and time consuming. The LANDSAT 
magnetic tapes provide continual coverage, but the user must have 
access to a computer (comparable to an IBM 360) and in most cases a 
computer technician, both of which are expensive. An Imager 100 system 
would cost as much as one million dollars. The only such system in 
Canada, is located at the Canada Centre for Remote Sensing in Ottawa. 

For these reasons it was impossible to incorporate any of 
those systems into the study. Consequently, available equipment had 
to be used (a density slicer and a color additive viewer). 

A major difficulty in mapping snow patterns from satellite 
imagery was snow drifting. This was especially true on the Cooking 
Lake Moraine. In its initial distribution, snowfall is more uniform 
than rainfall (McKay and Findlay, 1971, Western Snow Conference). 
"However, because it is easily moved by the wind, snow usually accumu- 
lates on the ground in a highly heterogenous manner, and it is, there- 
fore, difficult to obtain representative measurements" (McKay and 
Findlay, 1971, Western Snow Conference, p. 18). 

Physical features play an important role in snow concentration 
and distribution patterns on the Cooking Lake Moraine. Plage & isa 
photo of an open, flat, windswept section of Hastings Lake (Mt. Rose 
snow sampler indicating ll inches or 27.5 cm of snow). On the south 
side of Joseph Lake, trees provided a natural snow fence (56 inches 


or 130 cm of snow, Plate 3). During the winter of 1974-75 drifting 
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Plate 2 


Windswept Area of Hastings Lake 


Plate, 3 


Snow Drifting Patterns on the South 


Side of Joseph Lake 
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within and behind the natural snow fences was the bases for much of the 
Spring runoff. 
The snow pack tended to be fairly uniform in relationship 
to topographic and vegetation variation on the Cooking Moraine (author's 
field work during 1974-75). The amounts varied from north to south on 
the moraine, but the snow tended to be constant in heterogeneous zones. 
"The close relationship existing between snow cover depth and density 
and vegetation on the local scale persists even in regional patterns" 
(McKay and Findlay, 1971, Western Snow Conference, p. 19). McKay and 
Findlay further stated, "It is interesting to note that the standard 
deviations of the mean course value for each vegetation zone are small, 
considering the relatively few data used" (McKay and Findlay, 1971, 
Western Snow Conference, p. 19). . If this-is true, as it appears to 
be, it would mean a higher percentage level of predictability for 
snowmelt runoff on and around the Cooking Lake Moraine. 
The majority of the Cooking Lake Moraine is covered by a 

thick growth of poplars. 

Areas predominantly in aspen generally denote dry, well- 

drained sites. These may or may not be groundwater recharge 

areas. However, any snow accumulated in these areas has a 

high infiltration opportunity. Even under rather intense 

cultivation to remove this vegetation, the surface runoff 

from these areas is slight. If all were managed as recharge 

areas, the probability of improved ground water regime would 

be high'' (Swanson and Stevenson, 1971, Western Snow Conference, 

rope awe ts 
Swansen and Stevensons' report pertained to the eastern slopes of 
Alberta, a more humid region than the Cooking Lake Moraine. Yet, 


their suggestion to develop recharge areas to improve the groundwater 


regime could play a large role in an internal water management scheme 
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for the Cooking Lake Moraine, if depth of rooting was changed by 
changes in vegetative cover. 

Before an internal water management scheme could be developed 
for the Cooking Lake Moraine there are several important factors in 
need of consideration: Firstly, regional and local variations in 
annual precipitation patterns must be outlined using existing 
information services; secondly, background data on the effects of major 
storm tracks, topographical variations, and wind direction should be 
analyzed as to their effects on precipitation patterns. Chapter III 
provides a description of these factors as they appeared to relate 


to the Cooking Lake Moraine. 
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Chapter III 


REGIONAL AND LOCAL VARIATIONS IN TOTAL PRECIPITATION 


The emphasis in this chapter was placed on regional and 
local variations in total annual precipitation patterns. Maps, tables, 
and LANDSAT I imagery were used to identify the patterns as they 
exist on and around the Cooking Lake Moraine. 

The meteorological record for the Edmonton area is a 
moderately long one but it is for stations near, not within the moraine. 
A meteorological station was developed at Fort Edmonton, and obser- 
vations began on July 11, 1880 (Environment Canada, 1974). The aeerion 
was moved to several different locations and finally closed in 1942. 
After that time the meteorological data was obtained from the Edmonton 
Industrial Airport Station. The Edmonton Industrial Airport is on 
level ground surrounded by the city of Edmonton (Environment Canada, 
1974). 

Laycock summed up the meteorological data for the Edmonton 
station in the following; 

Let us review the 90-year record for Edmonton from 1883 to 1972. 

We find that the average water balance for areas with four inches 

storace (arain an average. value) is 17.6" Ppt. = (21.0% Pik. = 

4 Go ana. 2S. orby <0" Sts Change:)(Laycock,) 1973, app 2188). 
Laycock used Thornthwaite's water balance equation in this 


long term average water balance for the Edmonton station. With 
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Thornthwaite's water balance equation it was possible to calculate 
a balance among incoming precipitation, potential evapotranspiration, 
water deficit, water surplus, and soil moisture storage changes 
(Thornthwaite, 1955). 

A brief definition of the variables: 


Precipitation drizzle ~ rain, dew, snow, ‘elaze, and frost, 
which falls to or forms upon the earths surface. 


Potential Evapotranspiration = the total amount of water which 
might be transpired and evaporated by the vegetation with 

the energy supply available if moisture isn't limiting. 

Deficit = the difference between potential evapotranspiration 

and the amount of moisture that is actually available for 

evapotranspiration. 

Surplus = the amount of incoming moisture in excess of demand 
when soil moisture storage is fully recharged to root 
depth. This includes both surface and groundwater flow. 

Storage Change = plus or minus changes in soil moisture 

storage in the budget period. Changes in snow detention 

storage may be included. 

This formula was used to show water balance relationships for 
three soil moisture storage capacities (Table 8) on the Cooking Lake 
Moraine. The assigned soil moisture storage capacities used were: 

4 inch (10 cm) representing small grain crops and closely grazed 

pasture land, 6 inch (15 cm) representing mixed scrub forest and wooded 
pasture and the 8 inch (20 cm) representing mixed closed forest cover. 
The number of soil moisture storage capacity divisions will vary from one 
area to the next. For example, a 1/2 inch (1.3 cm) storage capacity 

may be used for intensively developed urban areas and 2 inches (5 cm) 


storage capacity for suburbs, where larger areas are devoted to lawns. 


A 10 inch (25 cm) capacity would be appropriate for more mature forest 
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cover. Almost all of the area under study can be classified using 
the 3 intermediate classes noted, because the urban areas are too small 
to be separately identified and dense, mature forest growth is present 
only in a few small groves. Refer to Thornthwaite, Laycock and Verma 
for factors involved in formulating soil moisture storage capacities 
for various soil complexes. 

To better envision climatic variations on the Cooking Lake 
Moraine it was beneficial to review the yearly water balance equations 
based on Thornthwaite procedures (1955, 1957) for five surrounding 
SGations (Plate 4, ‘Tables, 2,3, .4, 5; and 6)(metric conversions 
Appendix D). The stations were the Edmonton International Airport, 
the Edmonton Industrial Airport, Edmonton Namao, Camrose, and 
Vegreville CDA. The average soil moisture storage capacities, for 
the areas surrounding the stations were assumed to be four inches so 
that climatic differences might be stressed in comparison. The yearly 
water balance tables presented were computed for the four inch (10 cm) 
storage level. The records for the Elk Island Park station dated from 
1967 to 1975; but no year had a complete twelve month record. The 
data for the Ministik Lake station were also incomplete. Water balance 
tables for the Elk Island Park station were thus omitted. 

Because regional comparisons cannot be made on the bases of 
one station, there was little need to consider the meteorological 
data for the Edmonton station before 1950, except to show long term 
variations (Chapter II, Climatic Variations for the 90-year record). 
Several conclusions were reached based on analysis of the water balance 


equation tables for the five stations. The most apparent variation 
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Table 2 


Water Balance for Edmonton International Airport 1962-1974 


(Thornthwaite procedures using four inches storage)* 


Precipitation = UCPee =D) 1" Surpigs tstorage Change 
1962 Fa OR, aE aL ANE EN OSD Nea mera AY ~Le9 
1963 16.6 Se Ol eee veer): lO Csi hard 
1964 EBS) = (19.9 2-33.56 De -0 #33 
1965 2ae3 = (9.5 -20 .6:.9 & 4.3 = 49 
1966 20.2 25 AG wold oy anes rae: 0 ins peta =e 2 
1967 i oe ORC aa ate) at eS ial s) 
1968 ibe DAs, = Clo, r= 33'56- 24.0 =F), 
1969 13.2 =e CLO sa 0360.) 2 0 2 ut 
1970 Lien SGA a= eet tL. 9 ata ete 
1971 16.2 =O woh oo te 2. 0 = 5 
1972 Lao = (18.4 - 0 oa 2 4 agi ge 2 
eb ee: 23.0 = (G20 2) = 0 rere2 20 aah, 
Lore Lt iy ca Mat, U0 Re Bie Aes ie aaa 360 
Average 18.2 B se eh ye 1.6 2 


* See Appendices for metric conversions. 


Source: Atmospheric Environment Service Monthly Records. 
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Table 3 


Water Balance for Edmonton Industrial Airport 1950-1974 


(Thornthwaite procedures using four inches storage)* 


Precipitation 


Wa) 
20.4 
16.4 


15.8 
py 
7B ahs 
20.8 
17.8 


lI 


a a se? es 


(er 


(20.5 
O1gy,:5 
(ek s) 
C222 
hey) 
C229 
(22 
(2273 
(24,3 
(21.0 
C226 
C235 
C2252 
levy 
C2253 
2 
(21.8 
(COP as) 
a 
G22 eal 
(Zee 
(2 Decl 
Goi Hoses 
(Cs 


(PTO? 
C32", Se 53) 


- D) + Surplus *storage 
ei (pee: saad) er unee 
ani D aS ap aad bee, 
=o, 20) + lee, He ye 
Slt DP a eal Syed 
= 20) west if mb sd 
SURES ee IS oi Akrge' 
Soe a as fai. S 
1 Oius acre. O tno 
eee ots | ees - .7 
Sto. jar O cad gs 
pee ate O ats 
SL Ope tO - i.l 
4.6) 7+ 10 ay ee: 
9) tame = 8.'6 
me Oe): a 0 =. BAR) 
= e2 et 4.8 = 9 
AEG ce 0 ie oteid 
oP SIO) Bean ar og ews 
= 259) or. 0 - .9 
=i '5, 1) + 6 Ts Re: 
Oa leet - .4 
=O Smt he eh seo 
Pe Jet ee saee 
eG! )) fee 4) frites 
weed WAU Wh pene aleags) oe | 
me Oe nae 


ois, 


ee 


Atmospheric Environment Service Monthly Records. 
« See Appendices for Metric Conversions. 
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Tables4 
Water Balance for Edmonton Namao 1957-1974 


(Thornthwaite procedures using four inches storage)* 


Precipitation = (Pb. D) > Surplus Storage Change 
je 14.0 = G2R eS = 104") +0 a Si i 
L958 EPEg S22 = PF) 25 yt nats 
1959 Pout =e GS. 6 = Dra) i 1s - 4 
1960 liga key) SRC ee tm lce) att 6 o wih 64 
1961 Rib) Sue ese Mees a Pec aan | Ba eA 
£962 ig oe) See rs HO nO) var oo ae wares 
1963 1 eS: SF a PA ans Joes, rane 6 = .3 
1964 hsp) C2 = 7 ef reo Bid a Hem 
1965 Pas fees Bel.) 98 Or a, a - l.l 
1966 S72 ze CLOG aes Ol ect nies) ce) 
L907 14.0 =— G20R0F-) 8.9) a 5.5 ya hia 
1968 14.4 ee CON) oC) anes Gt =) ew 
1969 Mee, St 20S oer ae pel Me 
1970 i ee (2 = oot el. - 9 
1971 eye) = (P18 ea) oe 2 a es, 
L972 2062 =e AC LOCS a= Ee Orta 6 L =a! 
L973 22.6 = Vio = 2-66) <b 0 a Mbeya 
1974 Pte) =e VO20 one Ui) te ais - 1.8 
Average BS = C2 Ob = 2 ce Lady fea 


I 


Source: Atmospheric Environment Service Monthly Records, 


* See Appendices for Metric Conversions. 
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Sy, 
Tablen5 
Water Balance for Camrose 1951-1974 


(Thornthwaite procedures using four inches storage)* 


Precipitation a bees] Ot curpLus tStorage Change 
£O51 Noes cial leap! aca 0° by hale ela 2 PRA rhs) 
L952 3.) a ACME Un VaR) ent 8) A Nee a 
i me £O%.0 Wak (7 0 ia Say de  Bae wih 8 cine Be 
1954 22.0 ee (ee = OULi a Oo e4 a0) BY 
1955 16.2 Ae) su ey Foe i ee 
1956 16.7 ay AOL es cB hag Case =Hgiiaio 
L957 | esyeer =a Ouleat an eo) eae. t) see 
1958 11.4 —ec2lee <10. 740 + 0 
1959 16.5 2) AGO AO = sy) 2 6 eae 
1960 aie == -- 
1961 -- -- i 
1962 -- -- ae 
1963 -- ae -- 
1964 = = ce 
1965 ae Se = 
1966 Lore a= ACA ASmee SiS) am a8) z 6 
L967 eee as MOA ANS CRs) has 7 a 9 
1968 13.6 mal CAVA nies uses oy karen ol 
1969 £2. 4 ap gam Noreen aa a7 ce Ae 
Lice gl 16.9 = sa | Cop ae S) i 2.5 wil 2 
1971 14.3 SAL clean he Go - .4 
1972 20.6 = 420, 98= 352) 3.1 ay 
L373 28.9 aly iN ec Sar e 4M fee een es: is SI tae & 
1974 20:21 Son 10 ir a Siero 8 hey eae - 3.9 
Average 17.3 Sa ilies aes ee eed 
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Source: Atmospheric Environment Service Monthly Records. 


* See Appendices for Metric Conversion. 
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Table 6 


Water Balance for Vegreville CDA 1964-1974 


(Thornthwaite procedures using four inches storage)* 


Precipitation = CPL E.. = 'D)* + Strplts Storage Change 
1964 1320 wr 20 ele = Oo. <O Bg ee 
1965 ts. Ute tO) te = 5D 
1966 lee 2 —~ (20.4 = 8.9) 0 me 
L7o/ Lc Ue 20.20 9, Lb) tO coe 
1968 Tae 2 CLG Oey) ai ery 
1969 13.4 salem h G78 ew talent sent fir wal ope 
1270 Lath Fi WAPAVIAG IES SVE od 0 EL at's’ Cf 
BS ash Dae = "(2048 - 6.0) + 0 —aiw ears 
1972 ayers te LO eA ors )e te ed ee 
1973 758 Ws Se CADE ha NK, igadhen: 
1974 lesa s) Se Oso tei) pt 6G eras 
Average 1520 spel PAR CS alACes) Mae We eS) ti mae 


Sources: Atmospheric Environment Service Monthly Records. 


* See Appendices for Metric Conversions. 
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occurred at the Vegreville CDA station. Comparisons show that the 
Vegreville CDA station trailed the other stations in total precipi- 
tation by 2.3 to 3.2 inches (5.8 to 8 cm). Further comparisons showed 
that the Vegreville CDA station had the second largest deficit (-5.6 
inches or 14 cm). 

The largest average annual precipitation was received at 
the Edmonton International Mt BROLE (1825 inchesroni45.5 cm). The 
Edmonton International Airport station differed from the Vegreville 
CDA station in the total average deficit (3 inches or 7.5 cm smaller), 
surplus (1.5 inches or 3.8 cm larger), and storage change (.4 inches 
orl “em*larger). 

ParGlort thervartationgvin total tannual. prec ipitationswas 
caused by the location of the station, and by the type of gauging 
equipment used. The Edmonton International Airport, Edmonton Industrial 
Airport, and Edmonton Namao stations are first order stations (since 
1960) with nipher snow gauges. The Camrose and Vegreville CDA stations 
are second order stations and snow content is estimated to be one inch 
Ofawarer 10 tenmincnes Gf snow (2.5 em to 2) cm)... This could account 
for larger or smaller amounts of regional snowfall on and around the 
Cooking Lake Moraine than was indicated by average calculations for 
the surrounding stations. 

To compare meteorological data since 1961, for the five 
stations, it would be necessary to use the ten to one ratios for the 
snow data at the three airport stations. Due to the amounts recorded 


using a nipher gauge and the ten to one ratio, the actual precipitation 
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was probably greater than was indicated using the two gauging techniques. 

The Edmonton Industrial Airport station is located on level 
ground and is subject to a heat island effect by the city of Edmonton. 
This is indicated by the average potential evaporation calculations 
for, the station (22.3 inches or 55.8 em) which are higher than for any 
of the other stations. It was then assumed that variations in total 
precipitation would occur on the Cooking Lake Moraine. 

Three factors which may create variations in snowfall patterns 
on and around the Cooking Lake Moraine are: the paths of the major 
storms, topography, and wind direction. The first two factors, storms 
and topography, directly affect the amount of snowfall. The third 
factor, wind direction, tends to alter the distribution patterns of 


thessnow vy drifting it. 


Storms 


In July of 1964, Kozub, of the Department of Transport 
Meteorological Branch published a report entitled "Heavy Snowfalls at 
Edmonton", In this report Kozub classified 41 snow storms which 
occurred between 1944 and 1963 (Table 7). Kozub found that the 
majority of the storms fell into two categories: Surface Type I 
(12 storms) and Surface Type III (19 storms). The Gurface Type I 
storms were characterized by low-pressure areas near or on the British 
Columbia coast-line, which had been stationary or had moved in slowly 


from the Pacific Ocean (Figure V). The centers of this type of storm 
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4l 
Table <7 
Number of Snowstorms of Various Ranges, 
In Inches, By Months at Edmonton, Alberta, 


In the Winters 1944-45 to 1962-63 * 


Range Jan Feb Mar Apr May June July Aug Sept Oct Nov Dec Total 
Ze po eet 9 6 3 2 Z 2 5 8 Sy 
5.Ocm ot TEAS 

sheet Sipe) 9 a) 3 3 1 2 5 Z sl 
Teocm =) 9.6 

GeO = 49 2 2 2 3 72 3 3 4 24 
10cm -12.3 

5.0" (=. 5,49 5) 4 it 3 ia 
Oe.) eHiebic 9 2 1 ug 2 6 
Rocio 7c 

hed (MES eee Rea 1 2 1 4 
17.5cm-19.8 

S20 == 98.9 i 1 l 3 
20cm. -22.3 

9.0" = 9.9 l l i? 3 
22.5cm-24.8 

LO. Oa 1039 

Z5CM = 27.63 

11,0" -11.9 iL ih 
27.5cm-29.8 

1270" -12:59 

30.0cm-32.3 

£320" 2-135..9 

32.5cm-34.8 

bi ig ome Er 

35.0em-37.3 

15.0" =~ L359 

37.5cm-39.8 

E650" =-16.9 

40 -Ocm- 42 26s 

17.0" -17.9 1 1 
42.5cem-44.8 

18.0" -18.9 i 1 
45.0cm-47.3 

TOTAL ga Ati 7 Sh Ik an ed 2 3 (6 os A ca 


Source: Kozub "Heavy Snowfalls at Edmonton", 1964. 


* See Appendices for Metric Conversions. 
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tended to pass slightly south of Alberta. The Surface Type III storms 
were Pena ay eg by the movement of a low or a trough through 
Alberta at the surface (Figure VI).indicated by the arrow starting at 
Vancouver. With the maior snow storms, 31 out of 41, following 
these patterns it was apparent that these storms had a random effect 
On precipitation on and around the Cooking Lake Moraine. Kozub's 
study is only one of several acatrds with precipitation patterms for 
the Edmonton area. Refer to Longley, Whistance-Smith, Schrom, the 
Atlas of Climatic Maps, and Erxleben, and Potter for more in depth 
research pertaining to regional precipitation variations. It was 
impossible to depict any significant regional precipitation variations, 
for the Cooking Lake Moraine, from their research. Thus, the purpose 
of this study was to develop a more accurate means of estimating 
regional patterns of snowmelt runoff using remote sensing techniques. 
Effects of summer storms appeared to play a smaller role 
in precipitation variations on and around the Cooking Lake Moraine. 
Burrows, in a study entitled "Heavy Rainfalls at Edmonton", analyzed 
60 storms occuring in a period extending from April 21 - September 30, 
(1950 to 1965). “Surface, 850-mb, 700-mb, and 500-mb charts were 
carefully examined for each of the 60 storms of 0.75 inches or more, 
starting about 24 hours before commencement of rain at Edmonton" 
(Burrows, 1966, p. 5). Burrows found that an average of four storms 
a year with greater than .75 inches (1.9 cm) occurred during the 
summer season. Kozub limited the snowfalls he investigated to those 


depositing at least one inch on the ground (Kozub, 1964). Storms 
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Figure VI 


Paths of the Lows of Type 11 Storms 
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originating in the west produced the most frequent heavy rain storms, 
while heavy rain storms from the north occurred very seldom (Burrows, 
1966). Similar patterns were present in the major snow storm tracks. 
Burrows also found that the storms producing the most complete rain 
cover over a large area, developed from a cold low in the southwestern 
United States. Kozub did not indicate one type of storm which produced 
more complete snow cover than other storms. Refer to the previously 
indicated sources and to Reinalt and Longley for a more complete 


background. 


Topography 


Topography probably plays a smaller role in regional precipi- 
tation patterns on the Cooking Lake Moraine than storm tracks. Topo- 
graphy it appears, does play an important role in larger regional 
patterns. 

The Edmonton Industrial Airport lies 2200 feet above sea level, 
and the Vegreville CDA station is approximately 2084 feet above sea level, 
a very small variation. Hastings Lake, directly between the two stations, 
is 2418 feet above sea level, directly in the middle of the Cooking 
Lake Moraine (Figure VII). The Moraine peaks at approximately 2600 feet, 
20 miles (32 kilometres) to the south. Total precipitation patterns 
from the west to the east side of the Cooking Lake Moraine are probably 
affected by topographic variations. This may be partially 
illustrated by comparing the Vegreville CDA station's total annual 
precipitation record with the other four station's records. Other 


factors which affect regional variations ‘in precipitation patterns are 
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discussed in articles by Laycock and Longley. A method, which may be 
used to indicate the effects of regional variations in snow depth and 
water content patterns, on the Cooking Lake Moraine and surrounding 
area is presented in Chapter IV. 

Topographical variations play an important role in the 
duration of snowmelt runoff. The daily discharge data illustrate this 
in hydrograph form, and they also illustrate the significance of snow- 
melt runoff as compared to spring rains for the Whitemud Creek near 
Edmonton, and the Vermilion River near Vegreville (Figure VIII). 

The Whitemud Creek Basin encompasses 142 square miles 
(367,8 square kilometres) of bevelled till plain with lacustrine 
deposits in the north. The Vermilion River Basin (376 square miles 
or 973.8 square kilometres) is composed of much the same type of 
deposits with many more broad shallow depressions. 

It is apparent from the hydrograph that the Whitemud Creek 
has a much higner flow than that of the Vermilion River. The Whitemud 
Basin is a smaller system, which experienced a flashy runoff period 
after storage capacities had been recharged. The Vermilion River 
Basin experienced a similar peak flow, but due to the greater size 
of the basin and more depression storage the flow is extended over a 
longer period. The Cooking Lake Moraine (approximately 700 square 
miles or 1813 square kilometres) has approximately two or three times 


as much depression storage as the Vermilion River Basin (hummocky 
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moraine versus a ground moraine). Heavy tree cover and lack of a major 


creek or river system to drain the moraine impedes snowmelt runoff. 
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Figure VIII 


Hydrograph for the Whitemud Creek and the Vermilion River Basins 
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It is possible to estimate the eae runoff period as beginning 

a week later and continuing as much as a week longer than in either 

of the two smaller basins. Most of the snowmelt on the Cooking Lake 
Moraine, after soil moisture storage capacities are recharged, is left 


to be evaporated from ponds and land surfaces. 


The third factor which probably affects variations in total 
precipitation is wind direction. The prevailing monthly wind directions 
for January-April 1975, at the Edmonton International Airport, were: 
January West North West, February South, March East South East, and 
April South East (Environment Canada, 1975). 

Predominant winds, to some extent, do affect the localized 
snowfall concentration and distribution patterns on the Cooking Lake 
Moraine by induc ing drifting. Although, the majority of the drifting 
probably takes place from more localized winds. These winds are 
probably altered by topography and vegetative cover, thus causing a 
variety of snow drift patterns to occur on the moraine. This can be 
shown by reviewing a density sliced LANDSAT I Image. Density slicing 
accentuates pattern differentiation. patterns on the Cooking Lake Moraine 
(Plate 5. 

Major as well as minor drifting patterns can be detected on 
Plate 5. Note the following: the light gold areas, predominantly 
along the east-southeast borders of the lakes, represent areas of 


major build-up; the light blue on the lakes have shallow wind swept 
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Legend For Plate 5 


Localized Drifting Patterns Indicated on a Density Sliced 


Band 5 LANDSAT I Image of the Cooking Lake Moraine 


1 - cloud cover (blue) 
2 - minor snow build-ups (gold) 
3. - dense tree cover (black) 
4 - major snow build-ups (bright gold) 
5 - major snow build-ups at the forests edge (green) 
6 - very dense tree cover (red) 
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Plate 5 


a 
Naan, | Ca 


ae 


) ih 4) ante 
’ j arabe 
’ : ' th ; iN iit 4 { 
t S| 
A 
1 P ( 40 a. " 
”. \q it , 
, 
‘ | 
i _ ? = 
oN i ap 
1 if . h 
¢ va | ; ? 
y i : ; 7 
° a he { . + 
a 3 } : i ; 
4 i ; 
A " 
f i \ SS } : ( 
iu i i / 
ma eh vil ts as J 
i » 
fi + 
; |, ae ee 
, Caan ie 
i" i i 
i 4 ‘ j ‘ 
f j y bs 
5 iz | » 
ha : q i 
, i "EY Kok 
<=) 
‘ Sy 
) f 
‘ uy 
Ae ag a 
: » - 
‘ vere 
; : “ i ay 


| 0) 


onl 


snow cover; the medium gold density levels are minor build-ups of 
drifted snow; the black and red areas indicate thick tree cover (red 
the thickest); and the green, which borders much of the black areas, 
shows major snow build-ups on the edge of and in the first few yards 
of the forest cover. All of the above snowfall concentration and dis- 


tribution patterns are discussed in Chapter IV. 
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Chapter IV 


MAPPING OF SNOW CONCENTRATION AND DISTRIBUTION 


PATTERNS WITH LANDSAT I IMAGERY 


Four methods of Aer snow patterns with LANDSAT I 
imagery are discussed in this chapter. The first method involved 
the use of a black and white (scale of 1:1,000,000) band 7, 
LANDSAT I print combined with Levels I and II of Andersons' Land Use 
Classification System. The remaining three methods involved the use 
of Andersons' Level II categories, Thornthwaites' soil moisture 
storage levels, and a zoom transfer scope. The enhancement techniques 
used were: enlargement, density slicing, and tri-color projection. 
These procedures and conclusions based upon their use are discussed 
under five subheadings: 
Scale 1:1,000,000 LANDSAT I Imagery 
Scale 1:250,000 LANDSAT I Imagery 
Density Sliced LANDSAT I Imagery 
Tri-Color Projected LANDSAT I Imagery 
Water Balance Relationships and LANDSAT I Imagery 
The April 10, 1975, LANDSAT I imagery was not available, 
therefore, April 10, 1975, snow depth and water content data were used 
in conjunction with January 10, 1975, imagery. The reasons for this 


substitution are discussed in Chapter V under the subheading ''Problems 
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Encountered During the Study". 


Scale 1:1,000,000 LANDSAT I Imagery 


The standard form for LANDSAT I imagery is black and white 
positive prints at a scale of approximately 1:1,000,000. In this 
study black and white positive transparencies at an approximate scale 
of 1:1,000,000 were used for interpretation purposes. For illustrative 
purposes LANDSAT I Imagery, such as Plate 6, were presented in the form. 
of black and white prints which were made from the positive transparen- 
cies at same scale. The dark areas surrounding the lakes represent 
heavy deciduous tree cover on Plate 6. Lake boundaries were easily 
mapped, except on the west side of Beaverhill Lake and the southwest 
side cof Joseph Lake. In both instances the lake surfaces blended in 
with the surrounding snowcovered farm land. Tonal variations in the 
snow cover were almost non-existent at this scale. Light cloud cover 
was also apparent southwest of Cooking Lake. 

The land use patterns outlined on the first overlay (#1) 
were classified with Andersons' Level I Classification System (Anderson, 
1973). Only three categories were used, 02 - Agricultural land, 
04 - Forest Land, and 05 - water (Appendix F). The boundaries for the 
three categories were drawn without supplemental information (maps 
and low level photography), as Anderson suggested. Only general land 
use patterns for the Cooking Lake Moraine could be identified using 


these categories. Conceivably, Anderson's Level I Categories 03 - 
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Legend For Overlay 1 
Anderson's Level I Land Use Classification 
02 Agricultural Land 
O04 Forest Land 


05 Water 


Legend For Overlay 2 
Anderson's Level II Land Use Classification 
O01 Cropland and Pasture 
02 Lakes 


03 Mixed Forest Land 


Legend For Plate 6 


Andersons! Level I and Level II Land Use Classification 
on a January 10, 1975 LANDSAT I Image 
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Rangeland and 06 - Barren Land might have been used. These categories 
were not readily adaptable to the Cooking Lake Moraine for two reasons; 
firstly, there were only small scattered areas which represented the 
categories; and secondly, the sub-categories for each of the Level I 
categories were more representative of land use variations at lower 
latitudes. 

It was impossible to estimate runoff potential with Level I 
categories. It was possible, though, to make runoff potential estimates 
by combining actual ground measurements with classification based on 
Anderson's Level II categories by associating soil moisture storage 
Capacities with each category. The divisions used were: Agricultural 
land 02 - O1 = 4 inch (10 cm) capacity, Forest Land 04 - 03 (mixed 
scrub forest with small ponds and sloughs) = 6 inch (15 cm) capacity, 
Forest Land 04 - 03 (mixed mature stands with dense crown cover) = 
8 inch (20 cm) capacity, and Water 05 - 02 = total lake area. 

Once the soil moisture storage capacities and land use relation- 

ships had been estimated the next step was to calculate the total area 
encompassed by each storage capacity in the study area. The study area 
represents the section of the image that was density sliced for intensive 
study. The next step was to relate snow depth and water content data 
(taken with a Mt. Rose snow sampler) to each soil moisture storage level. 
This can be done for larger areas than the one represented by the study 
area, but once the area exceeds two square inches (12 square cm) dis- 
tortion on the image occurs. This is not a serious factor when the image 


is used only for land use classification purpose. If relationships, 


=) 
‘ rave i fi 


We wh ed 7 ii vd aie 


: lavas aoe oe. feta a: Alves gpa dom 03 


ae un 


os VS ae 0 Ren oe od uh aia Fy pein pau at 


sue bio acetal HE 


oe) 


= 
= 


ee aD ane ee ee, dons 
1 f DE te ' 


Ads otic) ae u pi: tial chine “th ne mi Nias wh Aan r - ‘a. 0'h 


Pe 


x 9% ute oiep Hae Oy a hid bey white SF don hae wt 


; 
hs Ga ae i rn y i s¥ ne % stl eore te pie Sah ut 


ay "i 
(av aah Ai my aan Cine? gw. te" yiedhes tia an em wie i . 


i A i ; 
i : f 4 ib oe 


Lol ie we er ty ae ho" ay 4 Lal fin | ‘vit owas ne. ae) i ; 
Fe he inh baht bin re su pena nanan ome Saat: aks 

lt sou u hit siadjotea oe rae vit ey nang ni 

6 ee va at i i aa wen ai at Supers + basis er a hieers aunt 
Sake hae es ie op | a ney 

b ese ea ur yin, ve i Bi ty 


: 
i 
i a: 


ThA 


si 4 eae oe ia oi, one rine oes tet 

aS veh} pial ie ag : is basi 4. we 

le, e age eae Perry boa Ned eth, 
. te 


Peal = = ht 4 ; Re ie | ite fae 7 
iy) 


aes i cosa ta l 
Dn ais ! 

vee 

ie 
‘aay 
ve 
¥ 4 bY f 
ie i i ” i 
ot See 
ie vod 

iat Ter, 

Pea Vt 


56 


though, are to be drawn between these images and density sliced images, 
to determine which areas are encompassed by each density level, serious 
errors in area percentages could occur. For discussion of image 
distortion refer to Photogrammetric Engineering and Remote Sensing, 
May edition 1975, for an article. by. Wong. 

Using this procedure, combined with Thornthwaite water 
balance tables it was possible to estimate runoff potential. These 
data were omitted from the discussion because more accurate calculations 
were obtained using a larger scaled image. This is discussed under 


the next subheading "Scale 1:250,000 LANDSAT I Imagery". 


Scale 1:250,000 LANDSAT I Imagery 


The January 10, 1975 LANDSAT I image was enlarged four 
times (sixteen times in area), from approximately 1:1,000,000 to 
1:250,000, to enhance the tonal variations in the study area. This 
was necessary to decrease error in area percentage calculations for 
each soil moisture storage capacity. At this scale the resolution 
was still good with only minor scan lines appearing on Cooking and 
Hastings Lakes (Plate 7). 

The enlargement was made from the same black and white 
positive transparency (approximate scale 1:1,000,000) as the print 
described in the previous section. It was necessary first to make a 
negative of the transparency so that a positive print could then be 


made. 


In this instance the sharper boundary distinctions made it 
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Legend For Plate 7 


Andersons' Level II Land Use Classification 
January 10, 1975 LANDSAT I Image 


Scale 1:250,000 


O01 Cropland and Pasture 
02 Lakes 


03 Mixed Forest Land 
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possible to determine finer boundary distinctions for each soil moisture 
storage capacity. Anderson's Level II categories were again used to 
relate soil moisture storage capacities to land use patterns. The 
divisions used were: Agricultural land (crop and pasture lands) 02-01 = 
4 inch (10 cm) capacity, Forest land (mixed scrub forest with small 
ponds and sloughs) 04-03 = 6 inch (15 cm) capacity, Forest land (mixed 
mature stands with a dense crown cover) 04-03 = 8 inch (20 cm) capacity, 
and Water 05-02 = total lake area. 

By using the procedures, for estimating runoff potential, 
outlined in the previous section, more reliable percentages were 
obtained with LANDSAT I imagery at an approximate scale of 1:250,000. 
Tables showing these results are presented and discussed under the 
subheading "Water Balance and LANDSAT I Trager! later in this chapter. 

One major problem that enlargement alone could not overcome 
was detection of the lack of uniformity in the snow cover. Consequently, 
it was impossible to detect the extent of snow drifting patterns in the 
study area or to make allowances for the drifting patterns in the 
Thornthwaite water balance tables. A method for overcoming this obstacle 


is discussed under the next subheading "Density Sliced LANDSAT I Imagery". 


Density Sliced LANDSAT I Imagery 


The density slicer used in this study was a Spectral Data 


Corporation model with a 32 false color level capability. This density 
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slicer is located at the Alberta Remote Sensing Centre, Edmonton, 
Alberta. The major advantages with the density slicer is opportunity 
to observe the image being split (step by step) into as many as 32 
false color levels. 

After testing several various false color combinations, of 
tnes3 2 density levels, the settings used were: ) 1-16, 28, and 29... These 
settings gave the best correlation between the density levels and the 
actual snow patterns on the Cooking Lake Moraine. These settings were 
employed throughout the study on all bands of LANDSAT I imagery. 

Where close inspection of standard LANDSAT I images did not 
show snow drifting patterns, density slicing did (Plates 8 and 9). For 
example the light blue area (#1) on Plate 8, band 5, represents snow 
which has drifted and has a high reflectivity (leaving a smooth, crys- 
talized surface) in the visible portion of the electromagnetic spectrum. 
The high reflectivity was illustrated by the cloud cover, in the lower 
left hand corner of the photograph, being recorded at the same density 
level as the snow cover. In contrast, Plate 9, band 7, showed several 
density variations within the same clouded area. Some of the variation 
was caused by photographic processing, rather than actual surface 
variations. Refer to the legend of each Plate for a description of the 
false color density levels, which are listed numerically on each 
photograph. 

To better envisage the density level variations the density 
slicer was adjusted to enlarge the section of the image representing the 


main study area (Plate 10 and 11). On Plate 10 for example, the various 
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Legend For Plate 8 


Regional Snow Pattern Variations Using a 
Density Sliced Band 5 LANDSAT I Image 


of January 10, 1975 


1 - Flat wind swept areas with shallow snow cover 

2 - Cloud cover 

3 - Areas where drifting is partially hampered by 
wind breaks 

4 - Accumulated snow cover from drifting 

5 - Forest edges and small openings where snow 
accumulated 


6 - Inner dense forest cover 
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Legend for Plate 9 


Regional Snow Pattern Variations Using a 
Density Sliced Band 7 LANDSAT I Image 


ofidanuany, LO oy > 


1 - Flat wind swept areas with shallow snow cover 

2 - Cloud cover 

3 - Drifting development induced by wind breaks 

4 - Accumulated snow cover from drifting 

5 - Forest edges and small openings where snow has 
accumulated 

6 - Inner forest cover 

7 - Haze or cloud cover of sufficient-density to 
impede reflectance 
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Legend For Overlay 1 


Soil Moisture Storage Capacity Levels 

4 - Four inch soil moisture storage capacity 
6 - Six inch soil moisture storage capacity 

8 - Eight inch soil moisture storage capacity 
(Based on text in subheading "Scale 1:250,000 


LANDSAT I Imagery'"') 


Legend For Plate 10 


Localized Snow Pattern Variations Using a 
Density Sliced Band 5 LANDSAT I Image 


of January) 0, 1975 


1 - Flat wind swept areas with shallow snow cover 

2 - Drifting partially hampered by wind breaks 

3 - Forest edges and small openings where snow has 
accumulated 

4 - Inner forest cover 

5 - Fields where drifting was at a minimum due to wind 
breaks 

6 - Accumulated snow from drifting 

7 - Haze or cloud cover existed of sufficient density 
to impede reflectance 
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Legend for Plate 1l 


Localized Snow Patterns Using a 
Density Sliced Band 7 LANDSAT I Image 


Of JanwaryLO,ge 7 


1 - Flat wind swept areas with shallow snow cover 

2 - Flat areas where drifting was hampered by wind breaks 

3 - Forest edges and small openings where snow has 
accumulated 

4 - Inner forest cover 

5 - Fields where drifting was at a minimum due to 
lack of wind breaks 

6 - Accumulated snow from drifting 

7 - Haze or cloud cover of sufficient density to impede 


reflectance 
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Plate 11 
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drifting patterns on Cooking Lake can be easily separated. The blue area 
represents shallow wind swept snow, while the medium gold and light gold 
colors represent areas where the snow has accumulated (light gold repre- 
sents the largest accumulations). The separation of drifting patterns 
were not as evident on Plate 11 as they were on Plate 10. For example, 
density level variations were more distinct on Plate 10 (band 5), than 
they were on Plate ll (band 7). This was especially true of the density 
levels on Cooking Lake. Further comparisons indicated that Plate 10, 
band 5, illustrated the best density distinctions, as they in turn 
related to the actual snow patterns cbserved during field work. One 
reason for the contrasts in Plate 10 and Plate 11 was the lack of open 
water or melting snow which would increase the tonal response in the 

near infrared (band 7). On January 10, 1975, the temperature was -24°F 
(-31°C) at the Edmonton International Airport, consequently, no melting 
was taking place. Dominant reflectance occurred in the visible portion 
Gf the Spectrum (band 5). “lt “shoulld be noted that on January 971975, 
the average wind speed was 15 miles per hour from the WNW at the Edmonton 
International Airport, causing fresh drifting in the study area. 

The approximate lake boundaries and snow sampling sites were 
outlined on Plates 10 and 11. Distinction of exact lake boundaries and 
snow sampling site locations was not possible because of two complications; 
the first complication resulted from the curvature of the TV monitor 
from which the 35 mm photographs were taken. The second complication 
resulted from the overlapping of similar density levels on the photo- 


graph, which in some cases made separate distinction of boundaries 
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impossible (for example, it would be difficult to distinguish a lake 
surface from a snow covered field if they were side by side). 

Using actual field measurements taken with a Mt. Rose snow 
sampler on April 10, 1975, it was possible to relate these readings to 
the false color density levels, Plate 10 was used as an example. 

The next step was to make an overlay from the January 10, 
1975, black and white LANDSAT I positive transparency (band 5), which 
could be overlaid on PlatelO (overlay 1). The overlay was constructed 
using a zoom transfer scope, which enabled the author to make the 
necessary boundary adjustments. The purpose of this overlay was to 
show the varying soil moisture capacities in relation to the seven 
density levels. A grid system was used to determine what areas were 
encompassed by each density level. The area distinctions were made to 
the Mearest wO,0f a square mile (1.6 square kilometres). 

These percentage measurements were not made with the instru- 
mentation available on the density slicer for two reasons: Firstly, 
the percentage readings represent the total area of that density being 
displayed on the TV monitor. The density levels cannot be broken 
down to show just the densities representing lake surfaces as Opposed 
to those representing farm land. Secondly, distortion occured along 
the borders of the image as it was displayed on the TV monitor. 
Distortion was caused by small cracks in the screening material, which 
allowed light to pass through and shadows to occur around the edges. 
Due to shadow interference on the sides of the photographs it was 


necessary to reduce the size of the study area from 350 square miles to 
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224 square miles (906.5 square kilometres to 580 square kilometres). 
The area is outlined and labeled "Study Area" on Plates VI and VII. 
Because each band of LANDSAT I imagery depicts a different 
portion of the electromagnetic spectrum (o4ute lal microns) slightly 
different tonal patterns appeared on each image. By overlaying and 
registering each band of imagery the complete tonal pattern for an 
area could be observed. A method of accomplishing this is discussed 


under the next subheading "Tri-Color Projected LANDSAT I Imagery'"'. 


Tri-Color Projected LANDSAT I Imagery 


"The Additive Color Viewer provides a simple and effective 
means Of analyzing and evaluating multispectral imagery" (12s, 1974, 
Data Sheet 9). The additive color viewer used in this study was a 
Spectral Data Corporation model (located at the Federal Forestry Building, 
Edmonton, Alberta). This was the only viewer available for use at the 
time this study was conducted. In the author's opinion the 12g Cor pora- 
tion's additive color viewer is superior to the Spectral Data viewer 
in every aspect and would have contributed more to the study. 

April 15, 1974, imagery in 70mm form, was used to illustrate 
a complication encountered with deep snow. The major complication was 
the lack of land use patterns that could be detected visually. Only 
band 4 and band 5 were available for the April 15, 1974 scene. 

On Plate 12 band 4 was falsely colored green and band 5 was 

falsely colored red; both were set on maximum light intensity. The 


fuzziness at the right and bottom of the photographs was a direct result 
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Regional Snow and Land Use Pattern Variations Using 


A Tri-Color Projected Image Combining Bands 4 and 5 


Band 4 false colored green - maximum light intensity 


Band 5 false colored red - maximum light intensity 
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Plates 12 
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Legend for Plate 13 


Regional Snow and Land Use Pattern Variations Using 


A Tri-Color Projected Image Combining Bands 4 and 5 


Band 4 false colored green - half light intensity 
Band 5 false colored red - maximum light intensity 
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of a lack of fine adjustment inet eunen cnt ten on the model of viewer 
used. 

Urban and forest patterns were enhanced with these settings 
(outlined on Plate 12 and Plate 13), while lake boundaries and tonal 
boundaries for snow cover were almost non-existent in the study area, 
making it impossible to extract pertinent information concerning snow 
patterns with those settings. The author found that snow exceeding 
30 inches tended to cause a washing out of tonal response in all forms 
of LANDSAT I imagery used except when enhanced by density slicing. 

Plate 13 represents the same bands as before with a different 
filter scheme. Band 4 was filtered green and set on half light intensity, 
while band 5 was filtered blue and set on maximum light intensity. Tonal 
boundaries were improved over those in Plate 12 (note forest patterns) 
allowing more information about land use patterns to be extracted. Even 
so, the usefulness of Plate 12 and Plate 13 in this snow study was 
questionable. 

The problem of deep snow camouflaging land use patterns could 
have been eliminated with an additive color viewer in which proper registry 
could be obtained. With proper registry, enlargements could have been 
obtained to increase the tonal variations in the land use patterns. Den- 
sity slicing of all of the bands used in the tri-color projections and 
selecting the band which best depicts the desired variations would be 
the next step (Plate 14). (Note the reduction in the number of density 
patterns.) If drifting patterns are constant from one year to the next 


the blue area would represent the heaviest snow cover. 
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Legend for Plate 14 


Regional Snow Pattern Variations Using 
A Density Sliced Band 5 LANDSAT I Image 


On April 2 gauko io 


1 - Deep accumulations of snow 
2 - Minor accumulations of snow 
3 - Windswept snow 


4 - Dense tree cover 
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Following the procedures outlined in the previous section, overlays 
for either land use patterns or density levels which represented snow 
cover could have been made using the zoom transfer scope. 

By using the viewer in this manner a more accurate relation- 
ship could be drawn among land use patterns, storage level patterns, 
and snow patterns. With the additive color viewer the user has the 
benefit of all the bands of one image laid one over another, producing 


better, tonaledistinctionss 


Water Balance Relationships and LANDSAT I Imagery 


Comparisons showing the applications of a LANDSAT I image 
in its standard form (black and white, 1:1,000,000 print) and an image 
enhanced through density slicing were used to determine water balance 
patterns on a section of the Cooking Lake Moraine. Black and white 
LANDSAT I imagery, in its standard form was the basis for all inter- 


pretations made in the study area. 


Black and White LANDSAT I Imagery in Standard Form 


Anderson's Level II Classification categories were used to 
relate land use patterns to soil moisture storage capacities for the 
study area in the following manner: Agricultural land 01 = 4 inch 
(10 cm) capacities, Forest land 03 (mixed scrub forest with smail ponds 
and sloughs) = 6 inch (15 cm) capacities, Forest land 03 (mixed mature 


stands with a dense crown cover) = 8 inch (20 cm) capacities, and 02 = 
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total lake area (Plate VI, overlay 1). Approximately four percent 

(an estimation based on fieldwork) of the study area would fall into the 
2 inch (5 cm) soil moisture storage capacity. It was impossible to 
detect areas of 2 inch (5 cm) storage level large enough (greater than 
-3 Of a square mile or .8 of a square kilometre, half cell size on the 
graph) for mapping from the imagery and consequently this category was 
omitted. Water balance tables (following Thornthwaite procedures) were 
prepared to establish surplus and deficiency patterns for each storage 
level (Table 9 or metric Table 9A). The potential soil moisture storage 
levels (with snowmelt) were below the 4 inch (10 cm), 6 inch (15 cm), 
and 8 inch (20 cm) storage capacity levels in April, 1975, where average 
snowfall was present. The dry summer and fall of 1974 caused a depletion 
of moisture levels in each storage category. 

Storage level types and the total lake area in the 224 square 
milesor 360 square kilometres study area are shown in Table 10. These 
are related to snowmelt water content, the percentage area encompassed 
within each of the storage capacities, and the potential surplus of 
snowmelt waters. The snowmelt water content figures used were the means 
of the (Appendix D) readings obtained at each field site that pertained 
to the given soil moisture storage capacity. The potential snowmelt 
water surplus calculations, in acre feet, were the sums of the moisture 
available in each soil moisture storage type at the end of October (1974) 
and the snowmelt water available in each type area in April (1975) minus 
the storage capacity. The figure obtained represents the available 


surplus for each storage capacity. Finally the potential surpluses of 
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Table 9 


Monthly Water Balance For Edmonton Internat ional Air port 


June to December 1974 


13 


June July Aug. Sept. Oct. Nov. Dec. Annual 
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I 5.04 534 Sed TES, 1.04 
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Def 
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Table 9A 
Monthly Water Balance For Edmonton International Airport 
January to June 1975 
Ree, Jan. Feb. Mar. Apr. May June 
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Table 9B 


Monthly Water Balance For Edmonton International Airport 
June to December 1974 


Deg bmee ace eal s Uk ve Pca a CCN ae NOV ee SP CG er. ina uMAES 
2G 14.5 Borel Pied 8./ Oa7 -2.8 -5.9 
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Table 9C 
Monthly Water Balance For Edmonton International Airport 
January to June 1975 (Metric) 
Jan. Feb. Mar. Apr. May June 
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snowmelt water for all storage capacities and the lake surfaces were 
calculated. It was assumed that soil moisture storage levels on the 
Cooking Lake Moraine would be equal to those at the Edmonton International 
Airport. 

The large variation in snowmelt water content from the six 
inch (15 cm) storage capacity to the eight inch (20 cm) storage 
capacity was a result of two variables: First there were many small 
openings (small pastures, small sloughs, and small ponds) combined with 
Sparse tree growth in the 6 inch (15 cm) storage capacity area, these 
combined to induce drifting in the 6 inch (15 cm) storage capacity. 
Secondly, the number of sampling sites for the 6 inch (15 cm) storage 
Capacity was very limited as compared to the number of sampling sites 
in the 8 inch (20 cm) capacity (refer to Appendix D). These calculations 
were also made assuming that average snow depth and water content measure- 
ments, for each storage level, would be representative of the area. This 
was not a totally correct assumption because the snow cover was not 
uniform. 

From Table 10 it was possible to calculate how much snowmelt 
water was contained in each soil moisture storage capacity on April 10, 
1975 and its relationship to the total area sampled. The objective was 
to determine what effects snow patterns had on snowmelt runoff. To do 
this the areas snow drifting patterns had to be mapped, and estimated. 
The potential yield of snowmelt water for the lakes was calculated 
without considering a different evaporation rate. Assuming this each 


lake would have risen five inches as a result of winter snowfall, aside 
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Table 10 
Potential Surplus Snowmelt Water Available in the Study Area 


INPADELLEOLELY 75 


October 

Snowme lt soil Potential 
Storage Water Moisture patent Surplus in 
Capacity Content Storage Lote Surplus Area Acre Feet 
qn ee Or aes son 34.3 7,785.6 
6" 7.6 - 46 8.06 2.06 Ss, £205.09 
oy Sie) 2.46 6.36 0 44.0 0 
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Table 10A 


Potential Surplus Snowmelt Water Available in the Study Area 


in April of 1975 Gnetrie) 


October 
Snowmelt Soil Potential 
Storage Water Moisture Of Suppiuspin 
Gapacity Content Storage Total Surplus Area Acre Feet 
10 cm ee Cm O cm [2537 Cm eee Cm 3433 = Non ea f 
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from runoff from land surfaces. This would have left lake levels 1.4 
inches (3.5 cm) below the levels at the end of April 1974, again without 
considering runoff from land surfaces. This statement was based on 
Thornthwaite calculations for potential evapotranspiration (18.4 inches 
or 46 cm) and precipitation (12.9 inches or 32.3 cm) for the period of 
May to October 1974. It should also be noted that there was probably 
a greater groundwater storage lag for 1974 than in 1975 because of 
increased precipitation during the winter of 1974. 

There would also have been greater amounts of storage recharge 
in the 8 inch (20 cm) storage capacity than the calculations indicate 
due to runoff from the 4 inch (10 cm) and the 6 inch (15 cm) storage 
capacity areas. A portion of the surplus snowmelt water available in the 
4 inch (10 cm) and the 6 inch (15 cm) storage aap iis shown in Table 10, 
would have been utilized to recharge the 8 inch (20 cm) storage capacity 


supplies. 


Density Slicing 


Under the subheading "Density Sliced LANDSAT I Imagery", 
procedures for enhancing images of snow concentration and distribution 
patterns were discussed. These procedures were used to show snow 
pattern variations for soils of different moisture storage capacities 
(Plate 10, Overlay 1). 

Table 11 was used to show the total areas, snow depths, snow- 


melt water contents, and the amount of snowmelt water available for soil 
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Level 
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Snow Depth, Water Content, and Potential Snowmelt Water 


Taplecu. 


Available for Soil Moisture Recharge and Surplus 


Total Area 


2 aesie) Mile 


Sore 
18.0 
78.0 
31.0 
32.6 

6.0 


in Each Color Density Level 


OnwADE LL EO Lo7S 


Snow Depth 


9.4 inches 


A) 
24.9 
17.3 
ayy 
Zn) 
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Table 1LA 


Water 
Content 


2.8 inches 


4.2 
7.6 
Shae 
4.1 


Potential 


78 


Snowmelt Water 


3,679.6 


7,526.4 
7,263.6 
16,307.2 
6,857.4 
Ike rally acs 
De, 99.95 


Snow Depth, Water Content, and Potential Snowmelt Water 


Available for Soil Moisture Recharge and Surplus 


Total Area 


ee 


emus) Sq. km. 


88.6 
46.6 
Zkdep 
82.4 
87.0 
ay Be 


Snow Depth 


MBA) 
29.8 
62.3 
43.3 
38.0 
63.8 
46.8 


cm 


in Each Color Density Lével 


onrApr ate lO spll97 SeGmetric) 


Water 
Content 


Potential 


Snowmelt Water 


453.9 
928.4 
895.9 
2.011 a5 
845.8 
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moisture recharge and surplus in each of the seven density level classes. 
The snow depth and water content measurements were obtained by relating 
density levels (by color) to actual snow measurements at snow survey 
Sites for April 10,2975 (Plate 10). Appendix’E was used’ to show 
density levels, snow depth and water content relationships. 

The next step was to relate the predetermined soil moisture 
storage capacities for each land use (taken from subheading "Scale 
1:250,000 LANDSAT I Imagery'') to the various density levels to determine 
percentages of areas in each. Percentages were obtained by using the 
Same grid system as before. Table 12 is a summary of these inter- 
relationships. 

Table 12A was used to show potential surplus snowmelt water 
for each storage capacity and the lake surfaces by image density levels. 
Surplus figures were obtained by adding the moisture available in each 
storage capacity at the end of October (1974) with the snowmelt water 
available in each density level and then by subtracting each storage 
capacity from each total. The figures obtained represent the available 
surplus for each storage capacity by density level. Again the final 
step was to calculate potential acre feet of surplus snowmelt water for 
the storage capacities and the lake surfaces by density level (Table 12B). 

From Table 12B it was possible to infer that the 4 inch (10 cm) 
storage capacity areas were affected noticeably more by snowdrifting 
than the other two storage level areas. This is shown by the lack of 
storage moisture at the end of October (1974) and the lack of surplus 


in April 1975, thus the potential surplus snowmelt water shown in Table 12B 
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Table 12 
Areal Relationships of Image Density Level Classes 


to Soil Moisture Storage and Lake Classes 


Density (hevelome sorelOrcm 6) "-orslovcem = 8 "sors20ecm™ Lakes 


1 LO hs | Loe 17.17 
2 ZN i -- 8) 20d 
3 -- 63.2 eee, 1.6 
4 15.4 36.8 (Bkee: 15.6 
5 Z1SES -- 6.8 10.9 
6 LOSS ~~ 6.8 18.8 
7 3.0 -- 0 6.3 
Table 12A 


Snowmelt Water Plus the End of October's Soil Moisture Storage 


in Excess of Storage Capacities for Each Image Density Level 


Density 

Levels “sae on = LOscm Gl Otten Once POV Ore 20.cne «lakes cm 
1 0 0 0 0 0 0 Zoo 7.0 
Z eee 3) 0) 0 0 0 Ged PO. 
§ bare, 9 2.06 me hs) 2.00 SPS. 7.6 LOa0 
4 0 0 0 0 0 GO Ore 9.8 
5 Ort a2) 0 0 0 aot 10.3 
6 4.2 1025 Zeoo 6.65 2.66 6.65 8.2 20.5 
7 CMe 


aa 0.66 1 565 0.66 | es ole 6.2 Loe. 
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Table 12B 
Potential Surplus for Areas in Each Image Density Level 


and Bach  Sterage Capacity =“April 10, 1975 


Acre Hectare Acre Hectare Acre Hectare Acre Heetare 
Density Feet Matres Feet Metres. Feet Metres 7 Reet Metres 
Levels Or ee ecm Oe One or clan oc. ero, (2 Orcm . “lakes: wom 
1 0 0 0 0 0 0 oa fo Oat 
Zz 164.7 2053 0 0 0 0 24956 ERGO 7 <i 
3 0 0 LOU. 7 94.0 C306. on 2 ae LOD) Whe) 
4 0 0 0 0 0 0 121315 149.55 
5 OF H7 Biol 0 0 0 0 891.4 110.0 
6 3354.6 Ay IE Tsk 0 0 05060) LUE 4) 307 De0e 319 66 
7 21063 33.4 0 0 0 0 Tih Oe ee. oi 
Total 3879 23 478.9 TOU 94.0 DO) el eee os DOGe ew 1b L90 26 


Table 13 
The Comparative Surpluses Indicated Using Standard and 


Density Slicing Based Calculations 


Storage Capacities Standard* Density Sliced** 
4 inch O acre feet 3,879.3 acre feet 
6 0 hOLLT 
8 0 2,347.7 
Lakes Te 12.0 FOGG. 7 

Table 13A 


Comparitive Amounts of Surplus Snowmelt Water Between Standard 
and Density Sliced Calculations Assuming No Precipitation 


Variation in the Standard Column 


Storage Capacities Standard Density Sliced. 30. 
10 cm O hectare metres 478.9 hectare metres 
LS 0 94.0 

20 0 28,9 8 


Lakes 953,66 1190.6 


Ss 


*Assume no precipitation distribution variation or drifting. 


““Assume drifting and precipitation distribution variation. 
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represents drifting. From Table 12B it was possible to determine in what 
storage capacities and in what density level areas snowfall had accumulated; 
assuming that established snowfall concentration and distribution patterns 
would not change dramatically from January 10, 1975 to April 10, 1975 

(snow depth and snowmelt water content measurements are considered to be 

the only variables). 

The level of accuracy was lower than the calculations presented 
show. The calculations were computed with a calculator and rounded off 
to the nearest tenth (for the values available). 

The differences between standard water balance calculations 
assuming no precipitation variations and density slicing based calcu- 
lations, were precipitation variations were considered, are shown in 
fable 13%. 

By referring back to Table 9 and the storage capacity values 
for April 1975, it was possible to determine that there was no surplus 
water available in any of the storage capacity areas. The acre feet 
of potential surplus snowmelt water shown in the standard column for 
the lakes were calculated by adding together the precipitation from 
November 1974, to April 1975, and then proceeding to calculate acre feet 
of surplus water (evaporation was not considered in these calculations). 

Larger variations could have been expected in normal (50"' or 
125 cm) or above normal snowfall seasons. 

The largest variation occurred in the 4 inch (10 cm) storage 
level, 3,879.3 acre feet (478.9 hectare metres). Table 14 was constructed 


to show the uniqueness of this remote sensing technique in detecting 
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differences in snow patterns (primarily drifting). 

By combining Table 11, Table 14, and Plate 10 the author was 
able to determine exact areas, within the 4 inch (10 cm) soil moisture 
storage capacity, where surpluses would occur. For instance, measure- 
able snowmelt runoff would occur in 57% of the 4 inch (10 cm) soil 
moisture storage level areas (primarily due to drifting patterns). 
Thus, with this remote sensing procedure it was possible to estimate 
areas having surplus (again possible runoff) not only for the entire 
area Of each soil moisture storage capacity, but also.«for each portion 
Gf these areas. (Table 14). 

The creas was taken one step further to calculate snow- 
melt runoff potential for two small areas, one surrounding Joseph Lake 
and one surrounding Hastings Lake (17.6 square miles or 45.6 square 
kilometres and 19.3 square miles or 50 square kilometres respectively). 
Maps were prepared showing the soil moisture storage capacities for 
the areas surrounding each lake (Figures IX and X).° Overlays were then 
prepared to show the seven density levels as they related to each soil 
moisture storage capacity and to the lake surfaces (as shown in Table 14). 
These overlays were constructed with a zoom transfer scope located at 
the Alberta Remote Sensing Center. By using the overlays in conjunction 
with Table 12B (showing the potential acre feet of surplus snowmelt 
water it was possible to determine snow patterns for each lake and its 
surrounding area. 

The next step was to calculate the area encompassed by each 


soil moisture storage capacity as well as each density level (following 
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Table 14 
Potential Acre Feet of Surplus Snewmelt Water in the Four 


Inch Soil Moisture Storage Capacity Areas on April 10, 1975 


Potential Acre Feet of 


Snowmelt Surplus 

Density Level % of Area Water Surplus Snowmelt Water 
Ll 2061 260 ines. 5 0 0 

2 20.4 4.2 we 164.7 

3 0 Te? 0 0 

4 15.4 B69 0 0 

5 2129 foil ol 89.7 

6 Loe Ce 4.2 3,354.6 

Z 3.0 Gee Die 270.3 

Table 14A 


Potential Hectare Metres of Surplus Snowmelt Water in the 10 cm 


Soil Moisture Storage Capacity Areas on April 10, 1975 (metric) 


Potential Hectare 


Snowmelt Metres of Surplus 
Density Level % Of Area Water Surplus Snowmeit Water 
1 2051 750eem 0 0 
2 20.1 L035 es) 20.3 
3 0 0 9 0 
4 15.4 9.8 0 0 
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Legend for Figure IX 


The Potential Soil Moisture Storage Levels for the Joseph Lake Area 


4 inch (10 cm) 


6 inch (15 cm) 


8 inch (20 ecm) 


The 


[@>) 
| 


I! 


Crop and Pasture Land 


Il 


Mixed Scrub Forest Land Combined with 


Small Sloughs and Ponds 


II 


Mixed Mature Forest with 


Dense Cover 


Legend for the Overlay 
Density Level Variations Representing Snow Depth and Water Content 


For the Joseph Lake Area on January 10, 1975 


Density level that represents 9.4 inches (23.5 cm) with 
2.8 inches (7 cm) water content 

Density level that represents 11.9 inches (29.8 cm) 
with 42) dnaches C10. Sicm) twater content 

Density level that represents 24.9 inches (62.3 cm) 
with 7.6 inches (19.0 cm) water content 

Density level that represents 17.3 inches (43.3 cm) 
with 3.9 inches (9.8 cm) water content 

Density level that represents 25.5 inches (63.8 cm) 
with 8.2 inches (20.5 cm) water content 

Density level that represents 18.7 inches (46.8 cm) 


with 6.2 inches (15.5 cm) water content 
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Legend for Figure X 
The Potential Soil Moisture Storage Levels for the 


Hastings Lake Area 


II 


inch (10 cm) = Crop and Pasture Land 


i 


inch (15 cm) Mixed Scrub Forest Land Combined 
with Small Sloughs and Ponds 


inch (20 cm) = Mixed Mature Forest with 


Dense Cover 


Legend For the Overlay 
The Density Level Variations Representing Snow Depth and Water 


Content for the Hastings Lake Area on January 10, 1975 


= Density level that represents 9.4 inches (23.5 cm) 
with 2.8 inches (7 cm) water content 

= Density level that represents 11.9 inches (29.8 cm) 
with 4.2 inches (10.5 em) water content 

= Density level that represents 24.9 inches (62.3 cm) 
with 7.6 inches (19.0 cm) water content 

= Density level that represents 17.3 inches (43.3 cm) 
with 3.9 inches (9.8 cm) water content 

= Density level that represents 25.5 inches (63.8 cm) 
with 6.2 inches /(20).5 em) water content 

= Density level that represents 18.7 inches (46.8 cm) 


with 6.2 inches (15.5 cm) water content 
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procedures outlined in the discussion for 12). Tables 15 and 16 show 
the results of these calculations. 

By combining the total calculations in Tables 15 and 16 with 
actual field measurements it was possible to calculate acre feet of 
surplus snowmelt water contained in the soil moisture storage capacity 
by density level. Tables 17 and 18 show these calculations for both 
Joseph and Hastings Lakes. 

Noteable variations in potential acre feet of surplus snowmelt 
water for the soil moisture storage capacities occurred between the 
Joseph Lake and the Hastings Lake areas. The largest variation in a soil 
moisture level occurred in the 6 inch (15 cm) capacity, 29.0 acre feet 
(3.6 hectare metres) for Joseph Lake and 148.4 acre feet (18.3 hectare 
metres) for Hastings Lake. Most of the variation was due to differences 
im land use patterns, which caused a variation in drifting patterns 
from Joseph Lake to Hastings Lake (farmland to forest land respectively). 
These variations had marked effects on snow drifting patterns in the two 
study areas. For example the trees surrounding Hastings Lake acted 
as natural snow fences, thus accumulating more snow than the 
agricultural land (4" or 10 cm) surrounding Joseph Lake. Yet, comparisons 
of the 4 inch (10 cm) at Joseph Lake to the 8 inch (20 cm) at Hastings 
Lake by density level, showed that more snowmelt surplus could be expected 
from the agricultural land (4.7 square miles) than from the forest land 
(7.4 square miles)(Tables 19 and 20). 

For Joseph Lake and the surrounding area, measurable surpluses 


occur in 3 density levels. For Hastings Lake and the surrounding area 
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Table 15 
The Percentage Area of Soil Moisture Storage Capacities 
in the Density Levels as Depicted on LANDSAT I Imagery 
for the Joseph Lake Area 


Density Levels A Orel em OiU0 Fag Cm Oh orezO.cm Lakes 
1 38% mile Ze ah 

Z 16 -- Ll 30 

| 5 28 14 bt 

4 34 (hE. 45 ie 

5 iis ae = ae 

6 il -- 9 20 

i -- -- -- iZ 

Table 16 


The Percentage Area of Soil Moisture Storage Capacities 
in the Density Levels as Depicted on LANDSAT I imagery 


for the Hastings Lake Area 


Density Levels PE or OF com OvnOnelucm | Shor 70 ci aka ee 
1 SibyA 16% 137%, 20% 

2 8 2 3 25 

3 del 21 17 l 

4 40 46 54 9 
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Taple Uy 


Potential Acre Feet and Hectare Metres of Surplus Snowmelt Water 


From Each Soil Moisture Storage Capacity Area Surrounding Joseph Lake 


Storage Levels Acre Feet Hectare Metres 
Ae Gr LO lcm LH. ine 
Dsor IL IeCm 2920 3.6 
olvor -20.cm 0 0 
Lakes LURES BE gta 243.3 
Table 18 


Potential Acre Feet and Hectare Metres of Surplus Snowmelt Water 
From Each Soil Moisture Storage Capacity Area Surrounding Hastings Lake 


Storage Level ent e Sah ACT eereer be Water eHectare Metres 
oor LO em 166.7 20.6 

Glory Loiem 148.4 io. 

hl ao eige7A ayes 1) 0 0 


Lakes 1,389.6 171.5 
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Table 19 


Potential Acre Feet of Surplus Snowmelt Water Contained in the Four Inch 


Soil Moisture Storage Capacity Area Surrounding Joseph Lake on April 10, 


1975 
Potential Acre Feet 

Density Level % Of Area Water Surplus of Surplus Snowmelt Water 
1 38 Vaso) 0 0 

2 16 4.2 «2 30.0 

3 5 7.6 3.6 168.9 

4 34 SiR, 0 0 

5 0 0 Ord 0 

6 ei 8.2 4.2 OM So yp 

7 0 0 0 0 

Table 1L9A 


Potential Hectare Metres of Surplus Snowmelt Water Contained in the 10 
Centimetre Soil Moisture Storage Capacity Area Surrounding Joseph Lake 


On Apr hi elO. blo 25 


Potential Hectare Metres 


Density Level % Of Aréale Water Surplus of Surplus Snowmelt Water 
1 38 7.0 cm 0 0 

2 16 10.5 ae SAG 

3 5 19.0 9.0 20.8 

4 34 9.8 0) 

5 0) 0) 0 0 

6 7 20.5 KO) 34.1 
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Potential Acre Feet of Surplus Snowmelt Water Contained in the Eight 
Inch Soil Moisture Storage Capacity Area Surrounding Hastings Lake 


Density Levels % of Area 


1 13 


onmApril 10,91975 


Water 

Zo 0 

Cen 0) 

ee 2-06 

Bie? 6) 

0 6) 

S.2 200 

0 0 
Table 20A 


Surplus 


Potential Acre Feet of 
Surplus Snowmelt Water 


Potential Acre Feet of Surplus Snowmelt Water Contained in the Twenty 


Centimetre Soil Moisture Storage Capacity Area for the Area Surrounding 
Hastings Lake on April 10, 1975 


Density Levels % of Area Water Surplus 
1 is PeOLCiMs 110 

2 5 10.5 0 

3 17 19.0 502 

4 54 9.8 0 

5 0 0 0 

6 13 20.5 6.65 

7 0 0 0 


Potential Hectare Metres 
of Surplus Snowmelt Water 


44.5 


43.9 
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measureable surplus occurred from only 2 density levels. The 8 inch 

(20 cm) storage capacity around Hastings Lake had only 59% (4.7 inches 
or 11.7 cm) as mich surplus moisture available as did the 4 inch (10 cm) 
storage capacity (7.9 inches or 19.7 cms) around Joseph Lake. This was 
primarily the result of land use variations (in drifting patterns) and 

a moderate snowfall (1975) combined with a dry summer and fall in 1974. 
Consequently, each storage level was kept below its full capacity. 

In the past it was neceSsary to have a net work of weather 
gauging stations, which would be representative of soil moisture storage 
Capacities. for the exact runofl potential calculations. “Yet) 1t was 
possible to show the variations in surplus on the entire moraine (Table 10) 
as well as in select areas (Tables 19 and 20) within the moraine. By 
combining remote sensing procedures, outlined -in this chapter, with 
Thornthwaite water balance tables (Table 9) it was possible to determine 
the effects of snow drifting patterns on potential surplus levels of 
the Cooking Lake Moraine (Table 12B). 

Under the present circumstances, with a substituted image, 
the water balance calculations presented are weakly based yet they 
represent a significant improvement over previous calculations made 
for the Cooking Lake Moraine. Other researchers who have used Thornthwaite 
water balance equations and tables to explain snow patterns in this 
region are Maciver, Landals, Wight, Erxleben, and Kakela. Other 
excellent sources dealing with runoff from snowmelt are the U.S. 


Army Corp of Engineers (1956) and Garstka (1964, refer to Bibliography). 
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Chapter V 


PROBLEMS ENCOUNTERED AND RECOMMENDATIONS 


The purpose of this chapter is to provide background information 
concerning the problems encountered during the study, to propose 
solutions to those problems, and to recommend several methods which 
would add significantly to further studies. The chapter is organized 
under two major subheadings: Problems Encountered During the Study 


and Recommendations for Further Snow Survey Studies. 


Problems Encountered During the Study 


The major problem encountered was attempting to obtain the April 10, 
1975 LANDSAT I image. On the morning of April 10, 1975 the author took 
snow depth and water content measurements on the Cooking Lake Moraine. 
It was a sunny, cloud free day. Accordingly, the author believed that 
a quality LANDSAT I image could be obtained. After sending eight 
letters (fotir to the EROS Data Center and four to the Canada Centre 
for Remote Sensing), making two long distance phone calls, and waiting 
seven months for replies, the author was informed that the image could 
not be obtained. Ms. Jean Heffernan, Manager of the Order Desk for 
Space Imagery, at the Canada Centre for Remote Sensing, called, saying 
that the April 10, 1975 LANDSAT I image was not available because it 


had been deemed of too poor quality to print and the magnetic tape had 
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been destroyed (Heffernan, phone call, 1975); with the magnetic tape 
destroyed, the image could not be obtained from the CCRS. 

Two days after Ms. Heffernan called, Mr. Leo A. Braconnier, 
Acting Chief, User Services at the EROS Data Centre, called, saying 
that EROS may no longer ship LANDSAT I or II imagery of Canada to 
anyone, anywhere without first ordering from Canada (Bracconier, phone 
call, 1975). Braconnier explained that Canada had sought to stop the 
U.S. from selling LANDSAT imagery of Canada (Brazil has sought the same 
type of injunction). Braconnier was not able to confirm the exact date 
that no additional imagery of Ganada could be sold, but did indicate 
that the April 10, 1975 images would not be shipped. 

The last option was to enlarge a micro-film image produced 
by Donald Fisher and Associates Ltd., Prince Albert, Saskatchewan. 
These images are not photographically corrected for distortion and are 
usually of poor quality. The image was "noise" (distortion) free 
ovér the study area, but_it was_of such poor tonal quality that the 
enlargement was not useable. Consequently, it was necessary to substitute 
the January 10, 1975 image for the Aprii 10, 1975 image. 

Problems encountered of lesser importance were: image quality 
deficiencies, inadequate tonal contrast on black and white prints, 
enlargement size, color brightness on 35 mm photographs, and the 
lack of gauging stations on the Cooking Lake Moraine. After acquiring 
the LANDSAT I imagery available for the study area (up to and including 
the January 10, 1975 imagery) it was found that EROS's and CCRS' image 


quality varied greatly. In every instance, the imagery from EROS was 
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of superior quality, clarity and definition. 

Tonal contrasts on a black and white print made from a 
positive transparency and color brightness on 35 mm photography created 
problems. An example of the problems encountered involved the December 
5, 1974 image used in Chapter I. It appeared from the print that the 
entire area was under a blanket of snow. In reality, there was less 
than an inch of snow on the study area (personal field work and snow 
pillow readings from the Edmonton International Airport). The technician, 
who printed the photograph, used too high a light intensity. The color 
prints were dull because of a low light intensity used during printing. 

The size of enlargements varied dramatically with the desired 
scale. A 1:500,000 print of an entire LANDSAT image is approximately 
24 inches (60 cm) by 24 inches (including the border), while a 1:250,000 
print is approximately 48 inches (120 cm) by 48 inches. In both cases 
the print would be too large to be included in the normal text of a 
paper. In snow survey studies where large areas are involved, overlays 
depicting land use and snow survey sites can be made from enlargements 
and photographically reduced to fit into the text. 

The lack of weather gauging stations (temperature and preci- 
pitation) on the Cooking Lake Moraine made it necessary to calculate 
Thornthwaite water balance tables with Edmonton International Airport 
data. This was an adequate arrangement, but not an ideal one. If time 
and funds did exist, a network of weather gauging stations would be 


beneficial in the study area. 
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Recommendations for Further Snow Survey Studies 


Although it is impossible for the average researcher to 
change governmental policies, there are several steps which can be 
taken. Firstly, a trin to the Canada Centre for Remote Sensing (or 
EROS Data Center) to familiarize oneself with the equipment, facilities, 
and key personnel who could aid in the study, would be beneficial. 
Secondly, a standing order for all magnetic tapes, as they become 
available, of the study area would save time and anxiety. This type 
of order can be cancelled for specific dates with a phone call if 
cloud cover was too extensive. 

The magnetic tapes should be incorporated into the snow 
survey. The tapes can provide a digital computer printout from which 
land use patterns can be related to soil moisture recharge levels. 
The key magnetic tapes (those just prior to the runoff) should be 
taken to the CCRS and displayed on the Imager 100, in a falsely colored 
digital format. The Imager 100 would enable a greater degree of 
accuracy to be obtained in the larger study areas. 

Intensive field work on the days of satellite pas sage is 
recommended. A network of predetermined snow survey sites should be 
set up by reviewing the previous winter's imagery with a density 
slicer. In this manner, snow concentration and distribution patterns 
could be established so that an appropriate number of field sites 
and locations could be determined. 

A skilled photo-technician, who has adequate laboratory 


facilities, is a necessity for any remote sensing study. The photo- 
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96 
technician can control definition, tone, and uniformity in all forms 
of photography. 

The final recommendation would be to have a low level flight 
(multi-band) flown just prior to the major snowmelt period. The flight 
should be aligned with the passage of the LANDSAT satellites. Complete 
coverage of a small study areais essential, in larger study areas, 
high level flights could be aligned with predetermined transects to 
incorporate as many snow survey sites as possible. This would give 
an excellent base from which interpretation could be made. 

If some or all of the previous recommendations would prove 
to be too costly, the procedures used in Chapter 4 of this thesis 
can be used. Complete preparation, execution of fieid work, and photo- 
graphic processing, as outlined in Chapter 4, would cost approximately 
$500 (excluding ealary and publication costs). Since the major 
objectives had to be narrowed from regional to local precipitation 
patterns it was only possible to indicate local drifting patterns 
and to estimate their effects on the water balance of the study area. 
If the study was enlarged to include regional variations in precipi- 
tation patterns the costs would increase proportionately with the 


size of the area being studied. 
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Chapter VI 


SUMMARY AND CONCLUS IONS 


Water quality and quantity problems are widely present 
on the Cooking Lake Moraine, making it unfit for many recreational 
activities. Several water importation schemes involving pipelines 
were proposed for the improvement of water supplies but these involved 
large cost for limited return. 

Alternative proposals involving better management of local 
water supplies appear to have promise, but they must be studied more 
thoroughly. Some of these involve management of snow distribution 
patterns using remote sensing techniques. 

Some of the factors which may affect snowfall concentration 
and distribution patterns on the Cooking Lake Moraine, were considered 
to be contributing factors in some snowfall variation patterns 
(relative to five local climate stations). The snow accumulation 
patterns from these storm tracks might be heavier on some portions of 
the moraine than other portions. Topographic differences may be 
responsible for variations in precipitation patterns, but without 
local long term weather gauging stations and appropriate imagery 
(late winter) it was impossible to establish patterns on a regional 
scale. Strong winds (especiaily WNW at the Edmonton International 
Airport) have contributed to variation in snow distribution in local 


as well as regional patterns of drifting. 
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Four basic techniques for determining snowfall distribution 
patterns were described. The techniques used were: Interpretation 
of scale 1:1,000,000 LANDSAT I Imagery, Enlargement to scale 
1:250,000 for LANDSAT I Imagery and more intensive interpretation, 

ensity licing of LANDSAT I Imagery for separation of snow depth 
categories, and tri-color projected LANDSAT I Imagery. The most 
accurate results were obtained by using a combination of LANDSAT I 
positive transparencies, density sliced images, field data, topo- 
graphical maps, and enhancement with a zoom transfer scope. The snow 
patterns (tonal) on LANDSAT I imagery at a scale of 1:1,000,000 were 
almost non-existant, but by combining enlarged LANDSAT I imagery 
with Andersons! Level II Classification System the potential for 
snowmelt runoff patterns was estimated. 

Enhancement by density slicing provided the best separation 
of snow concentration and distribution patterns on the LANDSAT I imagery. 
These snow patterns were mapped and related to potential snowmelt 
runoff variations for the study area, and for two more intensively 
studied parts of the study area. The runoff potential was also 
estimated for several soil moisture storage capacities. 

The major problem encountered during the study was the 
attempt to obtain the April 10, 1975 LANDSAT I images. The images 
were not available for two reasons: Firstly, the Canada Centre for 
Remote Sensing in Ottawa, had classified the images of the Cooking 
Lake Moraine for April 10, 1975 as being of too poor quality to print, 


consequently, they destroyed the magnetic tape. Secondly, the EROS 
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Data Genter, at Sioux Falls, Soern Dakota, may no longer sell LANDSAT 
imagery of Canada to anyone, anywhere. This arrangement was sought 
@atlier this year (1975) by Canadas ‘Substitution of the January LO, 
1975 images for the April 10, 1975 images was the best option available. 
Based on the problems encountered in this study three 
recommendations were made: Firstly, visiting the Canada Centre for 
Remote Sensing to familiarize the researcher with the equipment, 
Facilities, and key personnel that might aid in the researcher's 
study. Secondly, maintaining a standing order for all magnetic 
tapes of the study area, as they become available, could assure 
obtaining the image. This would speed up the processing of data 
for interpretation. Thirdly, having a low level flight (multi-band 
ideally side looking airborne radar) flown prior to the major snowmelt 
period (in conjunction with a LANDSAT passage) would ensure that useable 
imagery would be available. 
The purpose of this study was to develop a more accurate 
means of.estimating regional patterns of snowmelt runoff using remote 
sensing techniques. To fulfill the purpose of the study six objectives 


were formulated as the basis for research. 


Objectives 


The first objective was to establish what variations in 


winter snowfall are present in the Cooking Lake Moraine area. This 
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objective was not fulfilled completely for two reasons: Firstly, it was 
impossible to obtain the April 10, 1975 LANDSAT I imagery. This made it 
necessary to use January 10, 1975 LANDSAT I imagery in conjunction with 
April 10, 1975 snow depth and water content data. In making this 
substitution it was necessary to assume that snowfall concentration and 
distribution patterns would not vary appreciably from January 10, 1975 
topApril 10, 1975.5" Secondly,,snowfall ing l9/4-75.(to Jantiaryyl02n1975) 
was much lighter than normal and the major regional variation at that 
stage appeared to be due to drifting rather than to snowfall Suenencos: 
thus it was decided to focus upon drifting effects. With greater 
snowfall, more apparent regional differences might have been present 

and the procedure should not be discarded for use in other years. 

To develop a procedure using LANDSAT I imagery for deter- 
mining snow concentration and distribution patterns was the second 
objective. With the aid of a density slicer (using settings 1-16, 28, 
and 29), a zoom transfer scope, and positive transparencies of 
LANDSAT I imagery it was possible to fulfill this objective (for 
local patterns). A positive band 5 transparency was used as a basis 
for density slicing procedures and as a basis for the mapping of 
snow drifting patterns. The drifting patterns were mapped from 
density sliced prints with a zoom transfer scope. 

The drifting patterns were established and mapped on clear 
overlay material. The next step was to construct maps showing regional 
soil moisture storage capacity patterns. The overlays were placed 


over the (area surrounding Joseph and Hastings Lakes) maps and 
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relationships were established. Thus it was possible to determine the 
effects of snow drifting on a variety of scales with LANDSAT I imagery. 
These procedures should have been supplemented with more field data. 
The data should have been equally representative of all of the color 
density levels used. 

The third objective was to investigate the relationship 
between LANDSAT I images' he variations and snow concentration 
and distribution patterns. Comparisons of LANDSAT I imagery at 
scales of 1:1,000,000 and 1:250,000 showed little or no variations 
in) the snow cover. Thevmajor contribution of black and white 
LANDSAT prints was their use as base maps for determining land use 
patterns. Once the land use patterns had been determined, relation- 
ships were drawn between those patterns and the soil moisture 
storage capacities of the study area. It was then possible to estimate 
potential snowmelt runoff for each area on the moraine. 

To experiment with various enhancement devices to maximize 
discrimination of snow cover variations was the fourth objective. 
Three forms of enhancement were used in this study; enlargement, 
density slicing, and tri-color projection. Enlargement provided a 
better scale base map for distinguishing land use patterns on the 
moraine. Density slicing combined with snow survey data provided 
a bases for which snow drifting patterns could be mapped. Problems 
in registry rendered tri-color projections nearly useless in this 
study. If the registry could have been improved a useful image 


might have been obtained. 
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The fifth objective was to compare regional soil moisture 
recharge capacities and snow water content levels which were required 
to produce runoff. Long term water balance tables were calculated 
(Thornthwaite procedures) to show precipitation patterns for five 
weather gauging stations surrounding the moraine. They were: the 
Edmonton International Airport, the Edmonton Industrial Airport, 
Edmonton Namao, Camrose, and Vegreville CDA stations. Variations in 
gauging equipment (nipher gauges for the three Edmonton stations and 
stick measurements for the Camrose and Vegreville stations) made 
comparisons after 1960 difficult. 

Monthly water balance tables were calculated for the 
Edmonton International Airport station. These tables were used to 
determine the amount of water needed to recharge the soil moisture 
storage capacities. These data were compared to April 10, 1975 field 
measurements. Potential surplus water calculations were then made to 
determine the amounts of runoff which could be expected in each area. 
By combining the expected runoff data with maps derived from density 
sliced imagery it was possible to separate areas of different runoff 
potential from one another. It was also possibie to estimate potential 
surplus snowmelt water within a given storage capacity and break that 
data down into as many as seven different components (by color density 
level). This objective was fulfilled for localized areas. 

The final objective was to investigate the relaticnship of 
snow concentration and distribution patterns to soil moisture storage 


capacities so that runoff potential might be estimated. This was 
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accomplished for a 224 square mile (580 square kilometres) area of 
the Cooking Lake Moraine. By averaging the data obtained from the 
snow survey sites and relating that data, to the seven density 
levels, it was possible to estimate runoff. This procedure was 
also used on two small areas within the major study area, the area 
surrounding Joseph Lake and the area surrounding Hastings Lake. 

The major Ae eealey, en this portion of the study was 
the limited number of snow survey sites used, consequently, some 
density levels were represented by five or more sites, while some 
levels were represented by as few as three sites. Yet, with limited 
data it was possible to show runoff trends in the study area. 

Many snow surveys and investigations of drifting patterns can 
be completed utilizing the techniques and procedures described in this 
thesis. Particularly, it has been possible to show the effects of 
localized drifting in the Cooking Lake Moraine area. It should be 
possible to expand the mapping of snow patterns to a regional scale 
based upon the procedures outlined in Chapter IV. More complete snow- 
melt runoff relationships could also be determined. Future research 
should be undertaken to further refine the techniques and procedures 
upon which this thesis is based and to reinforce and expand its 


conclusions. 
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APPENDIX A 


Surplus and Deficiency Patterns for Edmonton 1883-1972, 
Westaskiwin 1921-1956 (33 years), and North Cooking Lake 1921-1929 (7 years) 
using Thornthwaite procedures (1948) with an average soil moisture 


storage capacity of four inches 


Edmonton ~- inches 
Surplus Deficiency Lake Level Change* 

A B C 
1883 4.0 O59 S049 oh RL F132 01 
1884 90 4.0 HO. 0 ri S600 - 6.0 
Le65. (On 4 Dao = OZ = 8759 = a0 
L686 0 De, = 8.9 a ee Be, ee) 
1507, (v0 8.0 =12.0 = 220) -12.0 
£o68, ~ 320 1.4 TAOS De: +1369 
£SOo). 0 14.9 -22.4 = 2294 a2 2at 
fegO* © 0 0 0 0 
LOO 1.0 Bete = B42 - 0.8 Tae 
Love, © 5.0 = feo = Al SO = 7.9 
1893 0.4 2.6 = Sto ren Peo ee ey: 
Looe ONS oe - 7.0 - 4,6 = eB S 2 
Povo. & Ose ok ~10.5 9e7 “49.3 
L670. Ose 7.0 = 9.7 + a733 = 4.9 
1897 9 8.5 =12.8 = 1255 -12.8 
LN ehe) <ammme xe, 10.7 =Lony - 14.5 = Loo 
1899 O 0 0 0 0 
1900 5.9 0 ss) S| 7 2050 +35.0 
1901 ips 0 TALS the ied amon a 
1902 “4.3 248 Th On THtBEO +2542 
P09.) ) (055 0 + 043 mi dsare 7 et 
1904 4.2 ee - 3.6 + aEOS.O +21.6 
0D, a0 6.3 - 9.5 =he OS =e0's5 
1906 0 7.4 = Vind =) SAL -ll.1 
1907. 4.0 0.7 2D 14.9 +26.9 
1908 §60 es ~10.8 10.8 -10.8 
1910 O Teo -11.0 - 11.0 -11.0 
roda-) 0 0 0 0 
Lote Onz - 0.3 - 0.3 - 0.3 
1913 Os 1.4 - 1.9 = 2363 - 0.7 
1914 0.8 0) ti Od rh died th De'O 
1915 0.6 as. ~ Sind - 1.4 + 0.4 
LOEGE abu 0.8 oe. 52 ot ies key et Lees 
1917 2.4 5.8 - 6.3 pane ae, Hob 
1916 0 3.2 - 4.8 - 4A.8 - 4.8 


hesesv <9 oseperabe et phe sie 


< ha Ma ity aay 


arn ‘ sae vine te sin 098) 


“ 7 
- TT nh biniahenas rites pom naan toni ah tee: tg ee 
tere tat a 
(j Pt ie 
"> \ ‘ Cut , 
Ys V.8 
ite .) eal i OL SF e 
GyElh eR 
i ey 
b. Soe be ad | ae: 
Pe Se ie N 


Soe a 
Bs. awl fat 7 

ae | so | toe x om 
Cfo ane 
Rue wl Tae 


+ 


Cat * ests 
8.51 | BSS Ee 
red Fits 


=)  ifhd 
vey page 


6 es whe 
oe hs SD Bike eae LL 


= ' 
kes he ” ; ee he + 
b Hav we 

e? va . t a) 5 de 
Var Ean 2S ee Mert ee Fe 
yok Sie «ane ek ce 
Shes” GO el ie ee ye Pe I 
. Ller OEE 

od ae 
Cates 


109 


(Appendix A cont inued) 
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(Appendix A continued) 


Surplus Deficiency Lake Level Change* 

A B C 
1a Lois eT OVE. Died -13.6 -13.0 -12.4 
£964 0 7ORS) -11.0 -11.0 -11.0 
L960 7 Ost 4.0 <2 59 +14.4 sa a0 lr/) 
1966 0 Joc - 7.8 - 7.8 - 7.8 
Evo) 0.7 7.6 -10.7 ~ 8.6 - 6.5 
B96e =O 6.1 Seley ey == eZ 
EMG. 0 4,7 - 7.1 sulel a trek 
LOFO.” O85 Shee) - 5.4 ee) - 2.4 
BOR O39 7.0 es OWS) =" 3.6 =1.0,..9 
Roya 325 ee oes 4 Igo ELS. 6 
EOS 


*Lake Level Change; A. Assuming no drainage from areas surrounding 
the lake; B. Assuming drainage from an area three times as large as the 
lake; and C. Assuming drainage from an area six times as large as the 
lake. Lake evaporation is assumed to be P.E. plus 1/2 D, thus the net 
level change without inflow would be 1 1/2 D-S. 


Source: Laycock, 1973. 
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(Appendix A continued) 


De ficiency Lake 
A 
14.8 =22 53 
bie) -1.8 
9.8 1 
20<0 ~30.8 
20 -30.0 
4.0 - 6.0 
18.8 =26..5 
8.3 ~ 9.8 
4.8 flo 
9.5 - 5.3 
1523 -13.0 
rigs) -11.3 
PS - 9.3 
16.0 =L5.o 
igs io - 8.8 
Tee} -18.8 
25 a AGS 
3.6 phd 46) 
9.0 -13.5 
16.5 -23.0 
p25 ~10.8 
9.0 - 9.8 
2 -20.0 
20.8 -31.3 
18.8 -28.0 
oss) ee 
js els. 0 
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APPENDIX B 


Available Reference Material Pertaining to 
Remote Sensing and Snow From the Canada 
Centre for Remote Sensing 


(Source: Canada Centre for Remote Sensing) 


Baker, D.R. Remote Sensing of Snow Fields from Earth Satellites, 
May, 1971. 


Barnes, J. et al. The Application of ERTS Imagery to Mapping Snow 
Cover in the Western United States, January, 1974. 


Barnes, J. Evaluate the Application of ERTS--A Data for Detecting 
and Mapping Snow Cover, May, 1973. 


Barnes, J.C. Satellite Surveillance of Mountain Snow in the Western 
United States, June, 1969. 


Barnes, J.C., Bowley, C.J. Use of ERTS Data for Mapping Snow Cover 
in the Western United States. March, 1973. 


Benson, C.S. Snow Cover Surveys in Alaska from ERTS-1 Data, March, 1973. 
Bryan, M.L. Fresh-Water In Interpretation from ERTS Imagery, 1973. 


Byrd, R.C. et al. Snow Measurement Using Millimetre Wavelengths, 
September, 1972. 


Carlson, R.F. Break-up Characteristics of the Chena River Watershed, 
Central Alaska. Dec., 1973. 


Carlson, R.F. Break-up Characteristics of Chena River Basin, July, 
1973s 


Deselm, H.R., Taylor, T.W. Vegetation Boundaries on ERIS-l Imagery. 


Dmitriev, A.Y. et al. Practical Use of Aircraft Gamma-Ray Survey of 
Snow Cover in the USSR. Sept., 1972. 
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Edgerton, A.T., Strogryn, A., Poe, G. Microwave Radiometric Investi- 
gation of Snowpacks. July 1971. 


Forestry Report, Canadian Forestry Service Northern Forest Research 
Centre, Remote Sensing, Earth Resources Technology Satellite. 
October, 1973. 


Golding, D.L. Snow Cover and Melting Snow from ERTS Imagery. June, 
LO7Ta. 


Grasty, R.L., Loijens, H.S., Ferguson, H.L. An Experimental Gamma- 
Ray Spectrometer Snow Survey Over Southern Ontario, Dec, 1973. 


Haefner, H. Status Report on the Current and Future Applications of 
R.S~ Technology in Switzerland, May, 1971. 


Hope, G. et al. Microwave Radiometry and It's Applications: Part 3, 
Hydrology, Oceanography. Feb, 1972. 


Howarth, P.J., Woo, M. A Preliminary Evaluation of ERTS-I Imagery 
for the Interpretation of Snow Gover S. Ontario, Aug., 1973. 


Itten, K. The Determinations of Snow-Lines From Weather Satellite 
Pictures. Sept., 1970. 


Kaminski, H. Snow and Ice in the N. Polar Regions Abolation 
Measurements from Sat. Photography, 1966-1972, March, 1973. 


Kennedy, J.M., Sakamoto, R.T., Passive Microwave Determinantions 
of Snow Wetness Factors. April, 1966. 


Kukla, G.J., Kukla, H.J. Increased Albedo in the Northern Hemisphere. 
February, 1974. 


Leeman, V. et al. The NASA Earth Resources Spectral Information 
System A Data Compilation. May, 1971. 


Linlor, W.I. Snowpack Water Content by Remote Sensing. Sept., 1972. 


Kondrat! Yev, K. Ya. et al. Spectral Albedo of Natural Grounds. 
suLy, “1906. 


Loijens, H.S., Grasty, R.L. Aiborne Measurement of Snow-Water 
Equivalent Using Natural Gamma Radiation Over S. Ontario. 
1973. 


Martiner, J. Evaluation of Air Photos for Snowmelt-Runoff Forecasts. 
Sepia Lr 
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Masry, S.E. The Analytical Plotter as a Stereo-Microdensitometer, 
ANG ALO Loe 


Meier, M.F., Applications of ERTS Imagery to Snow and Glacier 
Hydrology, May, 1973. 


Meier, M.F., Edgerton, A.T. Emission Characteristics of Snow and 
Ice in the Microwave Range, Dec., 1970. 


Meier, M.F. Measurements of Snow Cover Using Passive Microwave 
Radiation. isept., L072. 


Meier, M.F., Edgerton, A.T. Microwave Emissions from Snow--A Progress 
Report. May. (1971), 


McGinnis, D.O. Satellite Detection of Melting Snow and Ice by 
Simultaneous Visible and Near Intrared Measurement. Oct., 1972. 


Peck, E.L., Bissell, V.C. Evaluation of Snow Water Equivalent by 
Airborne Measurement of Passive Terrest-Gamma Radiation. Oct., 1971. 


Poulin, A.O. Hydrologic Characteristics of Snow-Covered Terrain From 
Thermal Infrared Imagery. Dec., 1973. 


Tavernier, R. Relationship Between Snow Cover and Discharge & Water 
Quality of Watersheds, Ardennes, Belgium, Nov., i973. 


Vickers, R.S., Rose, G.C. High Resolution Measurements of Snowpack 
Stratigraphy Using a Short Pulse Radar. Oct., 1972. 


Vonder Haar, T.H. Measurement of Albedo over Polar Snow and Ice 
Fields Using Nimbus-3 Satellite Data. Dec., 1973. 


Weller, G. Survey of the Seasonal Snow Cover in Alaska. July, 1973. 


Wendler, G., Carlson, R., Kane, D. Break-up Characteristics of the 
Chena River Watershed, Central Alaska. Dec., 1973. 


Wendler, G. Sea Ice Observation by Means of Satellite. March, 1973: 


Wiesnet, D.R. Detection of Snow Conditions in Mountainous Terrain. 
Sept. 1972. 


Wobber, F.J., Martin, K.R. Exploitation of ERTS-1 Imagery Utilizing 
Snow Enhancement Techniques. March, 1973. 
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APPENDIX C 


Taber ¢& 


Water Balance For Edmonton Namao A 1952-1974 


(Thornthwaite procedures using 10 cm storage level) 


Precipite tion 


LOOT 35.0 
1956 44.8 
|B Beye 45.3 
1960 47.3 
1961 44.0 
1962 MOSS) 
1963 28.8 
1964 hee 
1965 5L56 
1966 38.0 
1967 35.0 
1968 36.0 
1969 44.8 
19.70 43.8 
(Dow A he lars! 
Loy 2 sO ae 
Hes 56-5 
1974 54.5 

43.3 


Metric Conversions 


II 


(P.E.cm 


(53.8 - 
(55.3. 
(46.5 - 
(53.0 - 
(G33 -> 
(52. Se= 
(55.5 - 
(55.3.- 
(GiRere 
(49.0 - 
C52 ore 
(49.8 - 
(52.0 - 
(53.3 - 
(54. 5:00 
(49.5 


— 


- Dem) + 
25) + 
18) + 
Sup) ae 
Aes) ats 
Waa): 
16.5) % 
EW Eos Fate 
T9253) 1 
Ceo) as 
i EAS) We 
22S) ak 
es eat Wee 
Ly. 8) 
EC Bax 
ESSE) Bee» 
5) 35 


Surplus cm 
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Storage 
Change cm 


ara ery) 


aie me eas 


Source: Edmonton International Airport Monthly Meteorological Records. 
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APPENDIX C 


Table 2 
Water Balance For Edmonton International Airport 1962-1974 
(Thornthwaite procedures using 10 cm storage level) 


Metric Conversions 


Fs Storage 
Precipitation =  (P.E.cm - Dem) + Surplus cm Change cm 
1962 52.4 cm =) * 4058.0 2.4889 a 6.0 - 4.8 
1963 BL = (54.5 -14.8) er 1.9 = eo 
1964 44.0 = 1 (086. = 0208) 0 ER he | 
1965 DIVO ao (48 %80= 155) ou Oes tees: 
1966 40.8 a (47 =| 3) er a - 3 
1967 16 Pe) = C4975 = 193). eu ks, a 
1968 ayes) = (45.3 - 9) meen 8) einer! 
Pore 45.5 oe (48.5 - 9) rz 20 + 6 
1970 TAA) = GURS Sey) Ease ght amy FE) 
L971 40.5 2 Po Eells) Sect ae bh - 1.8 
1972 55 = (46 - 0) = ed 6 io3 
1973 SY eS ae (50.3 - 0) Lee) eos 
1974 46.8 = (5053) -"6.0) red Birr: - 9 
45.3 = CHO.G = J050R Fo 452 #38 


Source: Edmonton International Airport Monthly Meteorological Records. 
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1950 
1951 
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153 
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1956 
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L961 
L962 
1963 
1964 
1965 
1966 
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1968 
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Table 3 
Water Balance For Edmonton Industrial Airport 1950-1974 


(Thornthwaite procedures using 10 cm storage level) 


Metric Conversions 


*Storage 

Precipitation = (Po 6 .cm — Dem )a 4 Surp. us) cn Change cm 
323 GLH Bitte BB) scent 0 - a3 
54.0 aa (28.3°=' $5.5) + Bie: * 19,8 
41.0 a (5850) =) l5e0)ee “3 - 6 
63.6 5 (bua) 0) a yey S' 5 Ease 
49.8 = (49.8 - 0) cr To =m S 
DUeD = CS 7a eran oo 320 Gf ee 
DU. 3 ar (56. 08s) 16) ee 6.3 ogee 
3250 = (55.0 = 25) i 0 Ta para) 
143-54) a (60.8 - 18) a 220 - 1,6 
44.0 eg SDE a nO eae we eG) ee. o 
I) ae CaS gey et 5) as 0 aaa oft) 
5150 = (58.8 - 25) 1s 0 meee As. 
46.0 = (55s De melo) ara 259 - ps) 
Goes) se (59.3) =( 24.8) + ee mit ke) 
40.3 = (Oj Or— e200 junet 0 Tg Stes 
Soe, a (Sere lO a athe ee peer etre. 
38.8 = Lie 0 bo ye 0 - Pi 
3860 ra C56 .cn— oe ir Lea Loe 
eh ree, 2 C2505) = L deo | & 0 = 2a 
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(Table 3 Continued) 


tstorage 
Precipitation = (PsE.cm — Dem) + Surplus ‘cm Change cm 
LOW ately 5) = (593) =e 22) a Ie: at 8 
1972 50.8 = (53)583.=) 27) as 6.3 =wiroR 
Lo73 54.8 = (GAAS ae EES) Bs 0 Prd ge) 
1974 32.0 = Coe oe 6.3 - 3.8 
44.5 = (Gees Ae eas Dt ape Ab 


Source: Edmonton Industrial Airport Monthly Meteorological Records. 
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APPENDIX C 


Table 5 


Water Balance For Camrose 1951-1974 


(Thornthwaite procedures using 10 cm storage level) 


Metric Conversion 


Precipitation = + (P.E.cm - Dem) + Surplus cm 
LW 4243 = Coen 69S )) Das 
1952 3318 = Coe = Beco os 0 
£953 46.5 = (50.5886. 08 0 
B52 Bos) = (46.0 - 0) ii. Roe 
19155 40.5 = (51.8 - 13) + es) 
BRS: 41.8 a eR AR ee AC was L.5 
1937 37 08 = C2 Jie. 0) Gere 
1958 235) = (CsiAer yaa), a 0 
1959 Gi = G0s0 = Isso) er 10 
1960 -- = = -- -- 
1961) es =e ee oe 
1962 = = -- -- 
1963 == = -- -- -- 
1964 a Pte de ae =f 
1965 = = + -- ~~ 
1966 45.5 = (Sheets ae 0 
1967 S362 = Ss es Reni pte ie eG 
1968 34.0 = (50.5 - 14) ee) 


tstorage 
Change cm 


oe 


+. 


0 


hia 


ne 1) 


ede} Y, rn A nay 
ee Q 


let 4 


| hi ei i 


dhiss 


Nene 
a hrin | 


1969 


1970 


ESA 


WS) 722 


L973 


Lowe 


Source? 


Preci predation 


45.3 
aes: 
3550 
SASS 


(WA 


(Table 5 Continued) 


(P.E.cm - Dem) + 
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(54.0 ~ 14.8) + 
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Edmonton International Airport Monthly Meteorological Records. 
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APPENDIX 


Table 6 
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Water Balance for Vegreville CDA 1964-1974 


(Thornthwaite procedures using 10 cm storage level) 


Metric Conver 


Precipitation => ieGP. cm 
3455 =. 1 MSO: = 
46.3 = (50.3 - 
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Edmonton International Airport Monthly Meteorological Records. 
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The Number of Sampling Sites Used to Obtain Average Snow Depth 


APPENDIX D 


and Water Content Data For Each Density Level 
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Density Level Number of Samples Average Depth Average Water Content 
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APPENDIX: E 


Density Levels Ranked in Order of Snow Depth and Water Content 


Density Level Snow Depth Water Content 
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i 18.7 inches, 46.6 cm 6.2 inches; 15.5 cm 
2 Wo inche ss 420.oecm Ware ineiics, 10.9 cm 
5 iS.2 inches. Sevurem 4.1. inches, 20 .30em 
4 WP.3 Inches, 43.3 cm. 3.9 anchess 3.38 cm 
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APPENDIX F 


Land-Use Classification System for Use 


With Remote Sensor Data 


Levell 


Ol. 


02. 


03. 


O04. 


05. 


Urban and Built-up Land. 


Agricultural Land. 
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Forest Land. 
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Level. II 
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Residential. 

Commercial and services. 
Industrial. 

Extractive. 

Transportation, Communications, and 
ited Leateeie Soe 

Institut ional. 

Strip and Clustered Settlement. 
Mixed. 

Open and Other. 


Cropland and Pasture. 

Orchards, Groves, Bush Fruits, Vineyards, 
and Horticultural Areas. 

Feeding Operations. 

Other. 


Grass, 

Savannas (Palmetto Prairies). 
Chaparral. 

Desert Shrub. 


Deciduous. 
Evergreen (Coniferous and Others). 
Mixed. 


Streams and Waterways. 
Lakes. 
Reservoili 
Bays and 
Other. 
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06. Nonforested Wetland. 
Ol. Vegetated. 
O23 bare. 


074° -Barrenaband. 
Ol Sale flats. 
02. Beaches. 
03. Sand Other Than Beaches. 
O4. Bare Exposed Rock. 
05. Other. 


OS-- © liindra. 
OL ameltindras 


09. Permanent Snow and Icefields. 
OL. Permanent Snow and Icefields. 


—————_— 


Source: Geological Survey Circular 671. 
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